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ABSTRACT 

We present evidence of a >10-cr detection of the 10 /im silicate dust absorption feature in the 
spectrum of the gravitationally lensed quasar PKS 1830-211, produced by a foreground absorption 
system at redshift 0.886. We have examined more than 100 optical depth templates, derived from 
both observations of Galactic and extragalactic sources and laboratory measurements, in order 
to constrain the chemical structure of the silicate dust. We find that the best fit to the observed 
absorption profile is produced by laboratory crystalline olivine, with a corresponding peak optical 
depth of Tio = 0.27 ± 0.05. The fit is slightly improved upon by including small contributions 
from additional materials such as silica, enstatite, or serpentine, which suggests that the dust 
composition may consist of a blend of crystalline silicates. Combining templates for amorphous 
and crystalline silicates, we find that the fraction of crystalline silicates needs to be at least 95%. 
Given the rarity of extragalactic sources with such a high degree of silicate crystallinity, we also 
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explore the possibility that the observed spectral features are produced by amorphous silicates in 
combination with other molecular or atomic transitions, or by foreground source contamination. 
While we cannot rule out these latter possibilities, they lead to much poorer profile fits than 
for the crystalline olivine templates. If the presence of crystalline interstellar silicates in this 
distant galaxy is real, it would be highly unusual, given that the Milky Way interstellar matter 
contains essentially only amorphous silicates. It is possible that the z=0.886 absorber towards 
PKS 1830-211, well known for its high molecular content, has a unique star-forming environment 
that enables crystalline silicates to form and prevail. 

Subject headings: dust: extinction - galaxies: ISM - quasars: absorption 



INTRODUCTION 



Dust is a significant constituent of the observable universe, both impacting the appearance of galaxies 
and influencing their evolution. It affects the derived properties of some galaxies by attenuating the shorter- 
wavelength radiation emitted by stars and ionized gases, and re-emitting this radiation at longer wav elengths; 
more than 30% of the energy emitted as starlight is re-radiated in the IR by dust (IBernstein et al . 2002). 
Furthermore, dust influences physical processes in the interstellar medium (ISM) ranging from heating, 
cooling, and ionization, to those dictating the production of molecules and the formation of stars. Thus, 
understanding the nature of extragalactic interstellar dust grains is essential for characterizing the chemical 
evolution of galaxies, and for correcting the observations of local and high redshift objects used in studies 
spanning topics from galaxy morphology to cosmology. 

Despite the fundamental importance of dust in analyses characterizing galaxy evolution, however, there 
remain many open questions about the distribution and composition of dust in non-local galaxies, and about 
the possible cosmological evolution of dust properties. While significant evidence of dust-ens hrouded star- 



form a tion has been found in mid/far-IR and su b-mm studies of galaxies at moderate redshifts (jBarger et al 



1998 : Elbaz et al. IQQSl : Chary fc Elbaz 2001 ). detailed extinction curves have only been measured for a 



small number of local g roup galaxies, such as the Mfl ky Way (MW), SMC, and LMC, and for a subset of 



local starburst galaxies (Pei 1992; Calzetti et al. 2000l) . The local ISM can be studied both directly, through 



extinction and IR emission/absorption, and indirectly, via depletions inferred from gas-phase abundances. 
These studies reveal environmental variations in even local dust grain properties, such as differences in the 



— 

(Weltv et al. 


— „, , 
2001; 


Gordon et al. 


2OO3I 



3 Sofia et al. 2006^ . that are not fully understood. At higher redshift 



our understanding is even more incomplete, with limited direct knowledge about the detailed properties of 
dust in regular galaxies. The simplistic assumption that high-redshift dust grains are similar in size and 
composition to local dust may impact both derived galaxy properties, particularly in on-going large-scale 
galaxy surveys (e.g. SDSS, COSMOS), and techniques which are sensitive to the effects of dust atte nuation 
such as th e Type la superno va studies used to infer the acceleration of the cosmic expansion [e.g., Aguirrel 



(ll999a] lbD: lRiess et al.l (120041) ]. 



One technique which may be exploited to provide insight into the chemical composition of high-redshift 
ISM dust is the examination of absorption line systems in quasar spectra; particularly those in damped 
Lyman-alpha (DLA) absorption systems. These n eutral-hydrogen rich systems (log Nfff >20.3) con tain a 



substantial fraction of the neutral gas in galaxies (jWolfe et al.lll995t IStorrie-Lombardi fc Wolfd 12000), and 



are considered to be the best indicators of the chemical content of high redshift galaxies (jPettini et al. 



1994; 
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Prochaska et alJboOS : Kulkarni et al.lEoOSi) . Furthermore, since DLAs are selected solely by gas cross-section, 



rather t han by galaxy brightne ss, they pr ovide a relatively unbiased direct probe of dust in high redshift 
systems ( Ostriker et al.iil990: Fall et aljflQQS ) . Evidence for the presence of dust in DL A systems comes from 



Pei et al. 


1991; 


Pettini et al. 


1997) 



association of quasar reddening with absorptio n line strengths at l<z<2 (York et al. 20061). D ust in DLAs 



is also indirectly sugg ested by recent analyses (jKulkarni et al 



200E 



2007a, 



2011 



Peroux et al. 2006,, 2011i; 



Prochaska et al.ll2006f ) which find that the majority of DLAs exhibit low metallicities and lo w star-formation 
rates (SFRs). This f inding is in contradiction to chemical evolution model predictions ( Pei et al. 1999t 
Somerville et al.ll200l[ ) that are based on the cosmi c star- formation hist ory as inferred from galaxy imaging 
surveys (e.g. the Hubb le Deep Field. iMadau et all |l996h ). A scenario jPall et al.lll993l: iBoisse et al.lll998t 
Vladilo fc Perouxll2005[ ) in which the more dusty and more metal-rich DLAs provide greater obscuration of 
their background quasars, hence impeding their detection, could explain these findings. 

One of the primary open questions about high-redshift dust, which may be probed using DLA systems, 
pertains to the chemical composition and size distribution of the dust grains. Models for local dust generally 
assume a mixture of carbonaceous, silicate, and metallic/oxide grains with a given distributi on of grain sizes. 
In th e Milky Way, about 70% of th e core mass of interstellar dust is found in silicate grains ([Greenberg fc Li 
1999; Weingartner fc Drainel 2001 ). These silicate grains are produced only when oxygen is more abundant 
than carbon; the excess of silicate dust measured in the MW galactic ce nter has generally been attributed 
to a shortage of carbon-rich stars ( Roche fc Aitken 19851 Whittet 1987 ). At high redshift, however, most 
studies have focused only on characterizing carbonaceous dust, by probing either the 2175 A feature or the 
shape of the rest-frame UV extinction curve. Recent studies, by our group, examining quasar absorbers 
using the Spitzer InfraRed Spectrograph (IRS) though, have found clear detections of the 10 /xm silicate 



dust feature in these objects, with a >10-cr detection in the AO 0235+164 DLA system (jKulkarni et al 



2007bl I2OIII ) . These analyses suggest that the overall shape of the silicate absorption feature is best fit by 
an optical depth profile template derived from either laboratory amorphous olivine or from diffuse Galactic 
interstellar clouds, rather than from dense molecular clouds. These studies also find that the silicate feature 
is about 2-3 times as deep as expected from an extrapolation of the rio-Ay relationship for MW diffuse ISM 
clouds, suggesting that these systems may be probing dust in the inner regions of their respective absorbing 
host galaxies, rather than dust on the outskirts. 

In light of this recent success in detecting silicate dust in DLAs towards distant, moderately reddened 
quasars, we have now selected the 10 /im silicate feature in the PKS 1830-211 quasar absorption system for 
study. This particular quasar absorption system has been targeted for several reasons. First, the illuminating 
z=2.507 quasar is strongly lensed by a foreground (z=0.886), late-type (Sbc), low-inclination spiral, w hich is 
the presumed absorber host galaxy ( Winn et al. 2002; Courbin et al. 2002; Wiklind fc Combed 19981 ). This 
lensing produces two distinct images of the quasar; a brighter (NE) component with a magnification of 
5.9, and a second (SW) component which provides 1.5 times less magnification and exhibits signatures of 
dustier, more molecule-rich material (jWinn et al.l 12002). Towards the NE component and in the ab sorber 



disk Nhi ~2-3x102i cm'^, assuming T^pi„=100K (|Chengaluil Il999l : iKoopmans fc de Bruvnl 12005 ) wit h 
a density of n(iJ2) = 1700 — 2600 cm^"^ towards the more obscured SW component ( Henkel et al. 20091 ). 
While the spatial resolution of the Spitzer IRS is unable to separate these two lines of sight through the 
host galaxy, their integrated combination provides a significantly stronger quasar continuum detection, and 
resultant signal-to-noise r atio (SNR) for the silicate fea ture, than has been obtained for the previous analyses 
of silicate dust in DLAs (jKulkarni et al.ll2007bl I2OIII ). Second, while the majority of DLAs have relatively 
low amounts of dust, this source is part of a small subset of quasar absorbers found through radio/mm 
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surveys to be a molecule-rich, 21-cm (HI) absorber (|Chengaluij 119991 ) . indicative of an especially dusty 
and chemically more evolved quasar absorption system. Further evidence for the high dust content is 
provided by the relatively high measured absorber-frame extinction (Av<4 .1) and by the ratio of the total 



to selective extinction (Rv'=6.34±0.16) for this system iFalco et al.lll999l) . The high absorber rest-frame 



reddening of E(B-V)=0.57±0.13, with AE(B-V)=3.00±0.13 between the tw o quasar sight lines, is more 
typical of reddening values found for sub-m m sources [~0. 55, (^Calzetti*^2001 )1 than for e ither Lyman-break 



galaxie s [0-0.45 with a me dian of 0.16 at z~3: IShaplev et al. (2001); Papovich et al.l (|200ll )] or Mgll absorbers 
[0.002; York et al. ( 20061 )]. Thus, dusty absorbers such as this system may provide part of an evolutionary 



link in the sequence of SFRs, masses, metallicities, and dust content connecting the generally low-SFR, 
metal- and dust-poor general DLA population (detected in absorption) to the actively star-forming metal- 
and dust-rich Lyman-break galaxies (detected in emission). Third, the PKS 1830-211 system is known 
to exhibit mm/sub-mm absorption features revealing the presence of a multitude of atoms, such as [C 
I] iBottinelh et al.ll2009l) . and molecules, in many cases in multiple transitions, incl uding HCN, HCO"*", 
ortho-HaO, H NC, ortho-NH., N2H+ , H"CO+, H"CN, C S, CO, "CO, and C^^ O [e.g. IWiklind fc Comb^ 



( 1996L 119981) : 



Bottinelh et"al] J2OO9I): iHenkel et al.l |2008l ): iMenten et aP |2008l) : iMuller fc GuelinI |2008[ ): 



Muller et all poill )]. iMuller et al.l pood boill ) have found isotopic ratios {^^O/^^O, "0/i*^0, "N/i^ 



. , , . N, 

and ^^S/'^'*S) for this absorber which differ from those in the Milky Way. This abundance of molecules is 
suggestive of a chemically rich environment for the quasar absorption system, which we anticipate may be 
characteristic of systems which contribute significantly to the global metal content. 

We begin our paper with a description of our observations and data reductions in Sj2l and then in 
^ present our analysis utilizing the laboratory and observationally-derived silicate profile templates to 
estimate the optical depth in the 10 /im silicate absorption feature towards PKS 1830-211. We discuss 
in 21 tli6 implications of our finding that the prominently detected 10 /im silicate dust feature present in 
the PKS 1830-211 absorber has a crystalline silicate composition. Finally, we summarize our results in 
511 In Appendix El we discuss the impacts of variations in the adopted quasar-continuum normalization on 
our results, in Appendix |b1 we address two minor observational caveats which could potentially affect our 
analysis, and in AppendixOwe consider a broader range of observational and laboratory based optical depth 
profiles. 



DATA &: ANALYSIS 



2.1. Observations and Reductions 



Our observations, summ arized in Table [H were obtained using the InfraRed Spectrograph (IRS) on 
the Spitzer Space Telescope (| Werner et al.ll2004r ). These data consist of 69 individual spectral exposures, 
each observed in two nod positions, totaling 198 minutes (3.3 hours) of integration on-target, centered at 
(a, <^)j20oo=(18h33m39.89s, -21d03m40.4s). The two nods, obtained using the IRS Staring mode, facili- 
tate both background subtraction and the identification of spectral features induced by bad pixels. Given 
the spectral breadth of the silicate absorption features in previously studied quasar absorption systems 
( Kulkarni et al. 2007bL 2011 ). we chose to utilize the low- resolu tion mode (R'-^60- 130), rather than the high- 
resolution mode (R~600), in order to reach a higher sensitivity ( Houck et al.ll2004 ) . In addition to the LLl/2 
and SLl/2 spectra, bonus (third) order SL and LL spectra were obtained, which facilitate joining the first 
and second order spectra. Our complete data set spans a spectral region from 5.13 to 39.90 /xm; this broad 
wavelength coverage is crucial for the accurate determination of the underlying quasar continuum which we 



- 5 - 



require to measure the strength of the broad sihcate features. 



The IRS sht cannot individually distinguish the two lensed PKS 1830-211 quasar images; thus, we cannot 
probe two distinct lines of sight through the absorber. The lower-wavelength data (SL) are obtained with a 
3. "6-3. "7 X 57"slit and a l."8/pixel scale Si:As CCD, while t he longer- waveleng th data (LL) are obtained with 
a 10."5-10."7 X 168"slit and a 5."l/pixel scale Si:Sb CCD |Houck et al.ll2004l) . Given this resolution, we are 
unable to spatially differentiate the two lens ed images, which a re separated by (Aa, A(5)=(— 0."642, — 0."728) 
and embedded in a fainter Einstein ring (|Winn et al.l 120021 ). and instead we address, in this study, the 
composite spectrum which constrains the average dust prop erties within the quasar absorber. The NE 
lense d quasar image is significantly brighter at all wavelengths ( Lehai 200d Courbin et al. 2002 ; Winn et al 



20021 ). and we assume that our continuum is dominated by this line of sight. The weaker (SW) component 
is more heavily dust-ob scured, since it pa sses through a region of enhanced star formation in the western 
spiral arm of the host (jWinn et al.l 120021 ) , and will contribute more sign ificantly to o u r measured silicate 
dust absorption, as it does in detections of molecular absorption [see e.g. Ilehaij (l2OQ0l) : lFrve et al.l (|l997l) ]. 
Furthermore, we note that several neigh boring objects are spa tiall y inseparable from our quasar spectrum. 
These sources include [see illustrations in lCourbin et al.l (j2002r) and lWinn et al.l (|2002l )] an M4 dwarf star, a 
foreground galaxy at z^O.19, and possibly a second star and a second galaxy close to the z=0.886 absorber 
host. As explained in Appendix TO.2[ none of these other sources are expected to contribute significantly to 
our observed spectrum. 

Our analysis utilizes the Basic Calibrated Data (BCD) reduced spectra, and associated uncertainties 
and bad pixel masks, produced by version 18.7.0 of the Spitzer automatic pipeline. These spectra have been 
bias and dark subtracted, and a characteristic flat field has been applied. Additionally, the data have been 
corrected for systematic effects such as droop and dark-current drift, and cosmic ray hits have been identified 
and flagged in the corresponding masks. 

We have further processed the Spitzer BCD spectra by utilizing IRSCLEAN (vl. 9-2.0) to identify and 
remove rogue pixels and cosmic rays, and to reduce the background noise in the frames. In addition to 
replacing bad/rogue pixels identified by the masks and iterative algorithms, we have implemented a slightly 
non-traditional application of the IRSCLEAN software in which we also manually identified for replacement 
all pixels in the spectra backgrounds with visibly high/low signal in order to smooth the background. We did 
not apply IRSCLEAN -re^Aacevaenis, even for noted bad pixels, within the quasar spectrum or its environs, 
in order to prevent any spurious spectral signatures from being added to our data by the algorithm. We 
then manually removed bad/rogue pixels along or near the quasar spectrum using the epix task in IRAF, 
and replacing points with the median of the surrounding pixels. For pixels along the quasar spectrum, the 
unweighted average of bordering points along the spectral axis was used instead of the median. We were 
generally conservative with identifying and replacing these pixels, and only targeted strongly discrepant 
points, e.g. those pixels with saturated flux or isolated pixels with substantially higher/lower values than 
neighboring pixels. Pixels producing all other spectral features were retained as plausibly physical in ori- 
gin, and those which were produced by transitory events, such as cosmic rays, or pixel imperfections were 
effectively removed when we combined the data from the two nod positions. 

We combined the cleaned individual BCD exposures in each nod position using the IDL script coarf, 
which automatically propagates the uncertainties and mask files, and then subtracted the backgrounds using 
the IRAF task imarith, propagating the uncertainties in quadrature. For the SL data we subtracted one nod 
position from the other, in order to remove the background contribution. For the LL data, however, there 
exists an anomalous scattered light feature (see Appendix IB.1|) . This spatially-broad and intense feature 
prevented us from subtracting the two nods to remove the background. Therefore, we combined with coad 
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all of the ofF-target LLl frames (obtained while the quasar was in the LL2 order) and vice- versa to create 
off-target background frames, which were then subtracted from the spectra in each nod position. We applied 
the IRAF task fixpix to smooth over two weak foreground stellar spectra present in the ofF-target LLl spectra 
utilized for background subtraction. 

2.2. Spectral Extraction and Continuum Normalization 

We extracted the one-dimensional spectra for each order using the Spitzer IRS Custom Extractor 
[SPICE) algorithm, set for a point source extraction with optimal- weighting, using the default extrac- 
tion aperture in each order, with no subsequent corrections for fringing. Our final 1-dimensional spectra 
were extracted with a varying spatial width over the dispersion axis, as per the standard SPICE templates, 
sampling a broader spatial extent at the longest wavelengths within each order. A visual inspection of the 
two-dimensional spectra prior to extraction showed no apparent fringing in the SL2, SLl, or LL2 spectra, 
and only weak fringing in LLl. The IRSFRINCE manual stipulates that fringes are not spectrally resolved 
in the SL module, while in the LL module fringes resulting from one of the filters are present in the raw 
un-calibrated data, with an amplitude of 3-10%, depending on wavelength. These are largely corrected 
through the applied pipeline flat-fielding. As a test, we applied IRSFRINCE (v. 1.1) to our extracted LL 
data and found that the algorithm did not automatically detect any significant fringes, and when forced to 
correct for any measurable fringes, produced only a negligible statistical and visual improvement in the data 
quality, with small spurious spectral features resulting in a few cases. Therefore, we have applied no fringe 
corrections to our extracted spectra. 

In order to improve the signal-to-noise ratio (SNR) for our spectra, we have combined the one-dimensional 
extracted spectra from each of the two nod positions in an unweighted average. The fiux density uncertain- 
ties have been propagated in quadrature, and refiect both the measurement uncertainty estimated through 
the BCD pipeline reductions, and the scatter of the two nods relative to their mean. These uncertainties 
are dominated by variations between the two nod positions. Although the spectra extracted from the two 
nod positions are generally consistent, there exist regions with small discrepancies where, e.g., one of the 
nods exhibits a "convex" feature, while the other exhibits a "concave" feature. These small variations are at 
most 10% relative to the mean spectrum in the SL2, SLl, and LL2 data. In the LLl data, a discrepancy on 
the order of ±15-20% is observed over an extended spectral region from ^22-26 /im; this feature does not 
affect our measurements of the 10 and 18 /xm features in the z=0.886 absorber. After consultation with the 
Spitzer Science Center staff, we conclude that the origin of this feature is likely scattered light, as discussed 
in Appendix lB.il 

We have manually trimmed the nod-combined spectra to eliminate lower quality data points, and 
subsequently we normalized and joined the four SL-LL spectra to construct a single, continuous spectrum 
for the PKS 1830-211 system. These trimming and joining procedures are detailed in Appendix lA.ll Our 
combined spectrum has then been shifted to the rest frame of the quasar absorber host, assuming an absorber 
redshift of 0.886. The final combined spectrum is illustrated in the left panel of Figure[l] with the individual 
SL and LL orders differentiated by color. There, and in all other figures, we depict the absorber-rest-frame 
wavelength, unless otherwise noted. 

In order to search for and study the 10 and 18 /xm silicate features, we have normalized our spectrum 
by a fit to the quasar continuum. As detailed in Appendix IA.21 we explored a range of approximately 20 
different normalization fits, and identified a third order Chebyshev polynomial as producing the best fit to 
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the underlying continuum. This quality assessment was based on a combination of a visual inspection of the 
continuum fits and resultant continuum-normalized spectra, and an examination of the relative contributions 
by the continua to the reduced-chi-squared values for the template fits used in characterizing the optical 
depth and shape of the silicate absorption profile (for details, see iJ3]). Our final normalization is illustrated 
with a black line in the left panel of Figure [1] with the resultant normalized spectrum of PKS 1830-211 
depicted in the right panel. As is evident in this figure, there is a broad absorption feature present at 
^10 /im, as well as a shallower absorption feature near 18 fim. These features are associated with silicate 
dust in the z=0.886 foreground absorber, as discussed in the remainder of the paper. Additionally, there 
are several weaker absorption features in the 5-8 /zm region, and some emission features both longward and 
shortward of the 10 fj,m silicate absorption feature; the possible origins of these features are discussed in 
Appendix lA. 31 as they are not the primary focus of this analysis. 



RESULTS 



The z=0.886 absorber towards PKS 1830-211 produces prominent silicate absorption features near both 
10 and 18 fim, as is clearly illustrated in Figure [1] We have measured the equivalent width of the 10 ^m 
feature directly from thi s rest-frame, continuum-normalized spectrum, adopting a "nudge factor" of 0.3 
ISembach fc SavagelEooj) m estimating the continuum-fitting uncertainty. The total uncertainty associated 



with the calculated equivalent width includes the uncertainty from photon noise combined in quadrature with 
the uncertainty resulting from fitting the continuum. We find that the equivalent width of the 10 ^m feature 
is Wrest = 0.366 ± 0.029 ^m, resulting in a detection significance of 12.7-tT. If we exclude the contribution to 
the uncertainty from the continuum fitting, this significance would rise to a 28.4-cr detection. We have not 
estimated the equivalent width associated with the 18 /im absorption feature, since it lies at the edge of our 
spectral range and may not be fully covered by the current data set. 



3.1. Template Profile Fitting Procedures 

In order to understand the physical origin of the 10 /im silicate feature, and the chemical composition 
of the dusty material, we have fitted a series of template optical depth profiles (see Tables [2] and [3]) to our 
observed spectrum. Our fitting assumes simple radiative transfer through the cloud, such that I / Iq = exp[— r] 
where r = anTnorm', Tnorm IS the optical depth profile for the template, normalized to have a maximum 
peak depth of 1.0 over the full spectral extent of profile specified in Tables [2][3l The peak optical depth 
normalization factor (a„) for each template profile is selected to produce the minimum reduced chi-squared 
(Xr)- order to obtain the associated 1-a error bars, we have (cubic-spline) interpolated the two values 
in the Xr vs. r curve for which Xr ~ Xrmin = 1- The primary fitting has been performed over the 10 /im 
spectral region, i.e. 8.6-12.5 /im. This fitting region has been selected to focus on contributions to the Xr 
from portions of the spectrum near the absorption feature. For those objects in which the template profile 
additionally covers the 18 /im silicate feature, we have performed a second fit, termed the "full" fit, extending 
from 8.0-19.45 /im. This extended fitting region includes contributions to the Xr from not only the 10 and 
18 /im features, but also from the region of the continuum between these two features, which as discussed in 
Appendix [X] may also contain some emission from non-silicate sources. In those template profiles for which 
only a more limited wavelength range is covered, we have reduced this expanded fitting range appropriately, 
as detailed in Table IH 
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To constrain the composition of the material producing the 10 /im silicate absorption in the z=0 
absorber, we fit more than 100 miique optical depth profiles, derived from astrophysically observed and 
laboratory sources in the literature, as detailed in Tables [5J21 These templates probe a wide range of 
laboratory profiles for natural and synthetic terrestrial minerals, and environments including the Solar System 
(comets), circumstellar material, the Galactic diffuse interstellar medium (ISM), dense molecular clouds, and 
extragalactic objects (e.g. UltraLuminous InfraRed Galaxies; ULIRGs). For brevity, we include a subset 
of 18 of the most representative profiles in the bulk of this paper, as detailed in Table [2] (observational 
templates) and Table |3] (laboratory templates). However, in order to more fully explore the possible minerals 
comprising the detected silicate dust, we additionally examine an expanded subset of the templates, including 
44 laboratory and 18 observational profiles, including comets, in Appendix [Cl The full range of profiles are 
illustrated in Figure [5] (observed) and Figure [3] (laboratory). This complete set of templates spans a range of 
chemical composition s and degrees of silicate crystallinity, ranging from the purely amorph ous [diffuse ISM 
( Kemper et al.ll2004 . amorphous laboratory templa tes) to 10- 15% crystallinity [AGB star ([S ylvester et al 



20061: iKemper et al...2011i) ]. to 



IQQa iKemper et al.ll200ll: iMolster fc Kempeil 120051) . ULIRG (jSpoon et al 
purely crystalline laboratory silicates. 

The optical depth fits obtained using the 18 templates are presented in Table 01 We list the peak optical 
depth for fitting both solely over the 10 /im spectral region (8.6-12.5 /im) and over the 8.0-19.45 /im fitting 
range. In Figures S] and [S] we over-plot the fitted profiles on top of our PKS 1830-211 spectrum, using a 
template profile sampling of AA=0.01 /xm for uniformity. Considering our complete subset of 18 profiles, we 
find 0.11< Tio < 0.13 for the observationally-based profiles, and 0.10< no < 0.28 for the laboratory-based 
profiles. Most of these are visibly poor fits, as discussed in the following sections. Our best fit estimate of 
the peak optical depth is tio=0.27±0.05 for the crystalline olivine hortonolite template. This best fit was 
determined from both a visual inspection of the fits and based on the Xr i s-nd is at the highest end of the 
optical depth range derived using our template profiles. 

In order to ensure sampling uniformity for all of the computed fits, we have (cubic-spline) interpolated 
the values in each of the template profiles at the wavelengths at which our PKS 1830-211 spectrum has 
measured data points. By sampling every template at the wavelengths of our PKS 1830-211 data points, we 
ensure consistency when comparing the derived fits. This is essential since some of the laboratory profiles 
have been finely sampled with significantly sub-micron measurements, while the observationally derived 
profiles are more coarsely sampled. A consequence of sampling at our observed spectral dispersion scale is 
that we may potentially miss extremely narrow lines, but st udies examining un resolved atomic/molecular 
transitions in similar IRS spectra, e.g. star-forming galaxies (jSmith et al.ll2007l) , suggest that this is not a 
significant problem. 

Although the template profiles are sampled at a range of resolutions, since observationally-based tem- 
plates cannot achieve the high resolutions of laboratory measurements, we do not believe this is an issue 
for our analysis; the locations of the peaks play a more significant role than their shape. In order to test 
the effects of resolution variations, we have convolved three of the most-structured, crystalline olivine tem- 
plate profiles with two different Gaussians representative of the Spitzer IRS instrumental resolution; one of 
which represents the lowest achieved resolution typical of LLl, and one of which is more representative of 
an average spectral resolution across the wavelength range covering the 10 /xm feature. We find that for the 
more average spectral resolution there is a negligible change in the derived optical depth, with Ar=0.01. For 
the worst-case scenario, the optical depth is at most lower by Ar=0.04-0.05, which is consistent within our 
1-cr stated uncertainties. Since the instrumental resolution varies with wavelength over each IRS spectral 
module, it is difficult to correct for the instrumental profile. Given that our tests indicate no significant 
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impact on the derived optical depths, no significant improvements in the Xr ^oi our best-fitting profiles, and 
no visual improvement in the quality of the fit, we have made no corrections for the resolution variations 
between our template profiles and the PKS 1830-211 spectrum. 



3.2. Fits with Observational Templates 

While the laboratory-based profiles are arguably more precisely measured, and do not include the same 
observational uncertainties as astrophysical sources, we also include templates derived from astrophysical 
sources in this analysis for the following reasons: (i) they may contain compounds which are not commonly 
found terrestrially or synthesized in a laboratory; (ii) they may trace minerals in temperature and density 
environments more physically similar to our quasar absorber than can be achieved in ground-based labora- 
tories; and (iii) in order to establish the similarity or dissimilarity of higher-redshift extragalactic dust with 
that locally, it is essential to directly compare Galactic and other extragalactic dust sources with our sample. 



3.2.1. Fits with Galactic Source Templates 

We have first considered optical depth profiles derived from astrophysical sources within the Milky Way, 
including the diffuse ISM, dense molecular clouds, and circumstellar material. If extragalactic dust is similar 
to that in local systems, as is commonly implicitly assumed in extragalactic studies, then at least some of 
these sources should provide a good match to the 10 /zm (and 18 iim) silicate absorption in the PKS 1830- 



211 z=0.886 quasar absorber. Furthermore. iKulkarni et al.l (|201l[ ) have previously found that those quasar 



absorption systems in which the 10 /im silicate feature is not fit by laboratory amorphous olivine, are instead 
well-fit b y Galactic interstella r clouds, in which the dust is known to be composed of primarily amorphous 



material (jKemper et al.l 120041) . They have not found any instances wherein dense molecular clouds provide 
the best fit. 

In our analysis of the PKS 1830-211 absorber, however, none of these Galactic template profiles ade- 
quately fit the 10 /im feature. In general, these profiles exhibit three flaws: (i) the peak wavelength of the 
silicate absorption is shifted to shorter wavelengths relative to our absorber; (ii) they exhibit a single peak 
while we observe three distinct substructure peaks in the 10 /zm absorption region; and (iii) the breadth of 
our combined feature is relatively larger than can be reproduced by the Galactic templates. The best fit 
amongst the Galactic profiles is produced by a profile characteristic of an AGB star, perhaps the richest in 
crystalline silicates of our considered templates, but this is still a relatively poor fit when compared with the 
laboratory templates in the following sections. Neither the molecular clouds, nor the diffuse ISM profiles, 
nor the supergiant stellar environments produce a particularly good fit. It is remarkable that the Galactic 
molecular cloud templates do not fit the molecule-rich PKS 1830-211 absorption system; this may suggest 
that the molecular ISM in the PKS 1830-211 absorber differs from that in our own Galaxy. 



3.2.2. Fits with Extragalactic Source Templates 

We have also considered a silicate optical depth profile for an extragalactic source: a ULIRG (06301- 
7934, z=0.156), which may plausibly be more similar to our extragalactic z=0.886 absorber, despite the fact 
that it exhibits more robust star-formation signatures than are associated with our system. This ULIRG 
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was selected from a set of 12 such objects in Spoon et al. ( 20061) because it produce d the best match to 
our absorption system; coincidentally at 13% sihcate crystaUinity feem per et aL 2011 ) , it is also among the 
most crystalline-enriched ULIRGs from Spoon et al. ([2OO6 ) . 



We find that while the template profile based 
on this ULIRG produced a better fit than any of the Galactic sources, it still suffers from the same flaws 
in terms of an offset in the peak wavelength for the profile, and being unable to match either the observed 
substructure or the breadth of the feature. Removing the prominent 11.3 fim PAH complex emission feature 
from the ULIRG template, which is not seen in our system, improves upon the fit slightly but still does not 
allow us to reproduce our observed spectrum. 



3.3. Fits with Laboratory Templates 

3.3.1. Fits with Amorphous Silicate Templates 



We have next considered laboratory amorphous silicate profiles. As addressed i n iBowev fc Adamson 



(|2002bl ). the majority of analyses in both the local Galaxy, and at higher redshift [e.g.. ISpoon et al 



20061)], 

have considered the 10 ^m silicate feature to be produced primarily by an amorphous silicate, with residual 
structural f eatures attribut e d to ei t her c rystalline silicates or to non-silicate emission/absorption lines. The 
analyses bv lKulkarni et al. I (l2007bl I2OIII ) found that in most previously studied quasar absorption systems, 
the detected 10 /im silicate absorption can be reasonably reproduced by amorphous olivine. We have thus 
considered profiles representing both amorphous olivine {Al g2xF &2-2xSiO i) and pyroxene {M g^Fei-xSiO^) 
compositions, where < x < 10. We find that amorphous olivine provides a better fit than amorphous 
pyroxene, or than intermediate (not illustrated) amo rphous silicate profi l es, such as an amorphous version 
of the 21 Ferraro blend and an amorphous blend from Bowev fc Adamson ( 2002bl ). However, the amorphous 
silicate fits are still poorer than those produced by the crystalline silicates discussed below, or even by the 
more crystalline-rich observational templates. As with the observational template profiles, the laboratory 
amorphous silicate profiles are peaking at a lower wavelength, exhibit a shallower/narrower shape, and cannot 
reproduce the distinctive substructure we see in the PKS 1830-211 spectrum. If amorphous silicates are 
invoked to explain the 10 /zm feature in the system, we would also require significant superposed absorption 
and emission from other mechanisms. This possibility is addressed in fJH but leads to a poorer fit than our 
best case fit for a crystalline silicate template. 

We note that for both amorphous and crystalline (discussed below) olivine silicates, the values of r 
which are obtained fitting over the 10 ^m region are consistent with those obtained over the extended fitting 
range including the 18 ^m feature. This consistency suggests that our long-wavelength quasar continuum 
constraints are not unreasonable, although the relatively higher noise associated with the 18 /im data points, 
and the fact that we marginally cover the full extent of the 18 //m feature, places relatively weaker constraints 
on the derived optical depth measurements relative to the 10 /xm fitting. 



^Throughout our analysis, we note that amorphous ohvine/pyroxene refers to amorphous sihcates of an olivine/pyroxene 
chemical composition, since intrinsically all olivines and pyroxenes are crystalline in structure jHennindboTol'l . 
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3.3.2. Fits with Laboratory Crystalline Olivine Templates 

Given the relatively poor fits for both the amorphous silicates and for the observational sources, and 
the fact that the more crystalline-rich sources provide better fits than those which are crystalline-poor, we 
have next considered crystalline olivine silicates. We have examined fits for four such olivine templates 
{}Ag2xFe2-2xSi04) spanning a range of Mg/Fe ratios: synthetic forsterite (x=1.0), natural olivine (x=0.94), 
natural hortonolite (x=0.55), and synthetic fayalite (x=0.0). As illustrated in Figure El these minerals are 
able to explain the bulk of our profile shape and substructure, including the offset in the peak wavelength, the 
breadth of the feature, and the majority of the substructure peaks, without invoking secondary absorption 
and emission mechanisms such as would be required for amorphous silicates. The best fit is produced by 
hortonolite, both over the 10 /im fitting region, and when considering the expanded fitting region covering 
both the 10 /xm and 18 /xm absorption features, as illustrated in Figure [5] As is evident in Figure [Bl 
hortonolite well-reproduces the general shape of the profile, but cannot fully account for the leftmost peak 
(~9 fxm). Additionally, the rightmost peak (^11 /xm) is offset to slightly longer wavelengths relative to 
the PKS 1830-211 absorber spectrum. Comparatively forsterite, which produces the second-best fit, is able 
to match the shape and location of the rightmost (^11 /xm) substructure peak, but finds the middle peak 
(~10 /xm) and the intermediate ridge between these two features offset to slightly lower wavelengths. This 
strongly suggests that a crystalline olivine with an intermediate Mg/Fe ratio could simultaneously fit both 
peaks. While natural olivine is such an intermediate profile, it provides a poorer fit than either forsterite or 
hortonolite, fitting neither peak perfectly, suggesting an unexplored ratio between 0.55 and 1.0 may provide 
the best match. In order to determine the exact Mg/Fe ratio, longer wavelength data would be helpful as 
the similar p rofile shapes for th e four olivine compositions are significantly different longward of 20 /xm, as 
discussed by Jager et al. ( 1998h . 



3.3.3. Fits with Other Laboratory Templates 



Lastly, in order to explore options for a compound which could account for the leftmost (~9 ^m) 
substructure absorption feature, or which could perhaps explain all three substructure peaks simultane- 
ously, we have considered a range of other crystalline materials including an average serpentine, an average 
silica mixture, and two crystalline pyroxenes. We have also considered a SiC terrestrial/laboratory ma- 
terial, since although SiC has not been found in abundance in the ISM, it is found in meteorites and in 
evol ved stars and stellar p r oducts, including c a rbon stars, planetary nebu l ae, and some s uperri o ya remnants. 



Whittet et al 



(ll99ol):ISDeck et al.l Jl997l): Istanghellini et all feoOTh: Iziilstra et all |2006l ) : lYang et al 



|2004l) : lGruendl et all (|2008l ): iBernard-Salas et all |2009l ): lRho et all (|2010l) . We find that the best fitting of 



these alternative compositions is synthetic enstatite, a crystalline pyroxene, but that while it is an improve- 
ment over the other considered materials, and even over the laboratory amorphous silicate templates, it is 
not as good as any of the crystalline olivine templates. We also find that silica and serpentine, as illustrated 
in Figure [SI are successful at accounting for the leftmost (~9 /xm) substructure feature in our spectrum. 
Serpentine not only fits the '^9 /im feature, but also contributes some absorption in the 10 /xm region which 
is under-fit by some of the crystalline olivine profiles, such as the natural olivine and the fayalite. 
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3.3.4- Evidence for Multiple Minerals in Combination 



As our template fits have revealed that hortonolite produces the best fit, but that it cannot explain the 
•^9 /im absorption substructure, we have finally considered whether a combination of minerals could fully 
explain the absorption feature observed in the PKS 1830-211 spectrum. To test this we have considered 
a two-mineral mixture within the dusty quasar absorption system, combining hortonolite and a second 
mineral, such that I/Iq = cxp[— (a„ri^„orm + fenT2,norm)]- The optimal values for a„ and 6„ are obtained 
simultaneously, with the uncertainties conservatively determined by ascertaining parameter values which 
produce Xr ~ Xr mm ~ 1 within a given parameter, while holding the other parameter fixed. The resultant 
bi-variate fits are presented in Table [SJ and depicted in Figure [T) 



We find that the combination of hortonolite with average silica, average serpentine, or enstatite all reduce 
the Xr foi' the fit over the 10 /zm feature, and none of these profiles have significant features, for which would 
we need to account, in the 12.5-17 /^m region. We have isolated the two best combinations in Figure [51 in 
which it is evident that the strongest contribution is from hortonolite. The largest difference between these 
two fits is in the region between the 9 /im and 10 /im features where the enstatite combination exhibits 
an additional, very small, infiection. None of these fits quite match the height of the region intermediate 
between the 9 /im and 10 //m features, though, and while overall the fit is improved by the addition of 
the second mineral, there is only a weak contribution from this secondary material. For instance, in the 
combination of enstatite with hortonolite, enstatite is only contributing 10% of the column density relative 
to hortonolite. By adding these minerals we slightly reduce the hortonolite optical depth from 0.27±0.05 
to 0.24 at minimum (for enstatite). The optical depth of the second material is consistent with zero in 
every combination. We have also considered the combination of SiC with hortonolite, but find this is not an 
improvement. 

Alternatively, we note that it is possible that the '^9 /im absorption feature is produced by ammonia 
(NH3) absorption, w hich is seen in some ice- rich Galactic envir onments, su c h as W 33A, a dust-embedded 



young stehar object ( Gibb et al. 2000l 2004 ). As discussed in Gibb et al. (2004), NH3 has an umbreha 



mode transition at 9.3 /im which, when NH3 is a small component of the ice, shifts t o 9.0 /im. We have 
compared the modeled HsOiNHg (100:9) ice at 50K for W33A from iGibb et al.l |2000f ) with our spectrum 
and find that both the overall shape and location of that absorption feature is well-matched by our 9 /im 
absorption. Further more, ortho-NHs has previously been detected in the PKS 1830-211 z=0.886 quasar 



absorption system bv lMenten et al.l (|2008l ). with a total column density of 



In addition to considering combinations of hortonolite with a second material, we have also considered 
amorphous olivine in combination with either hortonolite or with SiG. We find that in the latter fit, SiC 
dominates the fit, and that while an improvement over pure amorphous olivine, is still a relatively poor fit 
in comparison with crystalline olivine, since the combined profile is unable to explain the shorter wavelength 
substructure features. The hortonolite and amorphous olivine fit is dominated completely by hortonolite, 
with little to no contribution from amorphous olivine. The amorphous olivine contribution produces a 
T <0.005, the lower limit considered in the combinat ion, leading us to c onclude that the contribution from 
crystalline silicates {\-^aml^cry) is >95%, following ISpoon et al.l (|2006l ). 



4. DISCUSSION & CONCLUSIONS 



In conclusion, based on our examination of a range of both observed and laboratory template profiles 
characterizing the 10 [xm silicate absorption feature in the obscuring z=0.886 absorber towards PKS 1830- 
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211, we have determined that our data are best fit by a laboratory crystalhne ohvine profile, possibly 
in combination with secondary silicate materials or ammonia. This finding is exceedingly unusual in the 
context of what is known about silicates in other astrophysical sources. The 10 /im and 18 /im features 
in Galactic interstellar matter are usually produced by the Si-0 stretching and 0-Si-O bending modes, 
respectively, of primarily amorphous silicate material. In the following section, we begin with a discussion of 
the plausibility of finding pure crystalline silicates in the PKS 1830-211 absorber, and then discuss several 
alternative scenarios which could potentially produce some of the observed structural features in the 10 /itm 
region. 



4.1. Crystalline Silicates 



Crystalline silicates have been observed in a multitude of Galact ic and extragalactic sources, but they 
generally contribute <15% of the silicate mass. iBowev fc Adamson |2002b) have examined whether the 
bro ad 10 ijui silicate feature c ould instead represent a superposition of numerous crystalline silicate features, 
but iMolster fc Kemper! (j2005l ) deem this scenario unlikely given that the corresponding expected crystalline 
silicate resonances at longer wavelengths have not been observed. We discuss in the following subsections a 
comparison of the PKS 1830-211 quasar absorber crystallinity with other Galactic and extragalactic sources, 
and with other quasar absorption systems, and finally we discuss the physical environments for crystalline 
silicates. 



Comparison With Other Astrophysical Sources 



While such high crystallinity as we tentatively observe towards PKS 1830-211 is very rare, objects with 
some degree of crystalline silicates are ubiquitous both in the local universe, and at higher redshifts. Galactic 
objects with noted crystalline silicates include comets (7-90% estimated crystallinity, depending on grain 
properties), interplanetary dust particles (IDPs), primitive meteoritic materials, young stellar environments 
including pre-main sequence stars such as Herbig Ae/Be stars and T-Tau stars, and evolved stars, particularly 
those with oxygen- and dust-rich outflows such as (high mass loss) AGB stars (10-15% crystallinity), post- 
AGB stars (20-60% crystallini ty in some post-AGB binary systems surrounded by a dust y circumbinary 



disc) , and planetary nebulae (IBouwman et al.. . 
200ll : lMeeus et allboOStlMolster fc Kempeilbooa 



Uchida et al 



2004 : 



20031: iGielen et al.ll201ll: iHonda et all 120031 : iKemper et al 



Wooden et al 



19991 ). Although a few 



rare stellar objects, such as the circu mstellar disk around the peculiar carbon star IRAS 09425-6040, exhibit 
up to 75% small-grain crystallinity (jMolster et al.ll200ir) . the material around the more populous low mass 
stars is generally less crystalline-enriched; the winds of low mass post-MS stars exhibit <10% crystallinity 
(jMolster fc Kempeiil2005l : iKemper et al.ll200lh. T h e diffuse ISM, at the lowest end of the spectrum, exhibits 



< 2-5% crystalline silicates (jKemper et al.l l2004t iLi et al 



20071). Among extragalactic sources, ULIRGs 



have been measured with at most 13% silicate crystallinity ( Spoon et al. 2006t Kemper et al. 2011 ). and the 
crystalline features are most pronounced longward of the 10 ^m region. These ULIRG templates provide 
a better match to our spectrum than any of the other observational templates. Galactic templates rich in 
crystalline silicates, such as the AGB OH/IR template, produce a better fit to our data than crystalline- 
poor templates, such as those representing the diffuse ISM. This all suggests that the PKS 1830-211 quasar 
absorption system at z=0.886 may indeed be unusually crystalline-rich. 
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4-1-2. Comparison With Other Quasar Absorbers 



Previous studies by our group of the silicate dust in other quasar absorbers did not find evidence for 
crystalhne structure, but no crystaUine-rich templates were utilized. Closer ex aminations of the fit residuals 
suggest additional unaccounted for substructure within the absorption features. iKulkarni et al.l (|2007bl . l2011l ) 
examined 5 quasar absorbers and found that the 10 /im feature was best reproduced by templates from either 
laboratory amorphous silicates of an olivine composition, or by the primarily amorphous diffuse Galactic 
ISM template. We have re-examined these fits in Aller et al. in preparation, using our expanded template 
library. We find that while in some of the systems, such as Q0235-I-164, the dominant component may be 
amorphous olivine, there are suggestions of additional crystalline silicate material in every system. In the 
case of the Q0937+5628 absorber we, in fact, find evidence that a high degree of silicate crystallinity may be 
present, similar to PKS 1830-211. Likewise, in 3C196 we find suggestions of significant crystallinity. These 
results lend further credence to our conclusion that the silicate dust in the PKS 1830-211 z=0.886 absorber 
may exhibit substantial crystallinity. 

Furthermore, it is plausible that the temperature of the absorber may impact the crystallinity. Measure- 
ments of the spin temperatures in some of these system s which exhibit the weakest sig natures of crystallinity, 
e.g. Q 0852-F3435 [T^ <536 K; ISrianand et all i2008h ] and Q0235-hl64 [T^=210 K; iKanekar fc Chengalur 
( 2003hl. may be lower than that within the PKS 1830-211 system [T^ ~1000 K in some of the H I material 
Subrahmanvan et al.l (1992)]. The primarily amorphous silicates in these other quasar absorption systems 
could be linked to their lower spin temperatures. (We note, however, that the spin temperatures in quasar 
absorbers are somewhat uncertai n, as there is debate a bout the fraction of th e quasar radio emission covered 
by the foreground absorber, e.g., Curran et al. ( 2005 ): Kanekar et al. ( 20091 )). 



4-1-3- Physical Conditions 



If the features in our spectrum are indeed produced by crystalline silicates, it would be interesting to 
understand both the stimulus and physical environment producing these crystallines and whether they are 
similar to the crystalline silicates detected in more local systems. Co nsidering the relatively large value 
of the total to selective extinction [Ry=6.34±0.16; IPalco et al.l |l999t )] estimated for the PKS 1830-211 



absorber, which is suggestive of large grains associated with relatively dense regions, it is possible that 
the environment in the galaxy may be conducive to grain growth which could be li nked with the apparent 
build-up of crystalline silicates detected. One should, however, note cautions by iMcGough et al.l (|2005l ) 
about the interpretation of Ky values in such lensed systems. In light of the measurements of crystalline 
silicates in other astrophysical sources, and in the laboratory, there persist a number of puzzling questions 
for the z=:0.886 absorption towards PKS 1830-211, if crystalline silicates are invoked to explain the spectral 
structures. 

First, it is unclear why there is no evidence of amorphous silicates in the 10 /im region of the PKS 1830- 
211 absorption system. Galactic sources, even those noted to be relatively crystalline-rich, are spectrally- 
dominated by the co ntribution from the amo rphous silicates in the 10 //m region, as is wel l-illustrated by the 
collage in Figure 4 of i Bouwman et aL ( 2001 ) and in the analysis of Bouwman et al. ( 2003 ). In our spectrum, 
however, we see significant potential crystalline structure, which could imply a deficit of the amorphous 
silicate material. It is not clear what mechanism would produce such a surfeit of crystalline material in the 
z=0.886 absorber towards PKS 1830-211, or alternatively, such a suppression of the amorphous silicates in 
the 10 /xm region. In the 18 fim region, where crystalline silicates might be expected to also be prominently 
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detected, our data are consistent with a broad amorphous-type structure with no distinctive crystalhne 
features; however, we note that our visibly noisier data in this region could possibly obscure narrow crystalline 
signatures. 

The second question which remains unclear is the temperature of the media producing and hosting the 
crystalline silicates. Galactic crystalline silicates are found in a range of different temperature environments 
and may be formed through di ffering physical proc esses. Temperature estimates within the z=0.886 absorber 
vary significantly according to iHenkel et al.l (j2008( ) , who estimate kinetic temperatures of ~^80 K for 80-90% 
for the system, with warmer temperatures for the remaining 10-20% of the material, some of which is ^600 K 
and generally concentrated in a spiral arm of the absorber galaxy. Detailed studies of laboratory crystallines, 
in combination with observations of crystallines in Galactic sources, indicate that efficient production of 
crystallines requires temperatures in excess of ~1000 K, and that if the material cools too rapidly after 
formation then the material becomes quickly amorphized. However, since we are observing the silicate profile 
in absorption rather than emission, this suggests that the material may be relatively cool, although we note 
that, in principle, some warm material could also be detected in absorption. The detection of crystalline 
silicates in a cooler medium is not completely unusual; for instance meteorites exhibit significant quantities 
of crystalline silicates. The presence of such cool crystalline silicates in an extragalactic environment would 
suggest that the crystalline material may have been transported in, e.g., a cooling outfiow from an initially 
hotter medium; in ULIRGs, the cooling outfiows stemming from evolved stars within the star-bursts are 
invoked to explain the temperature discrepancy between the observed crystalline silicate absorption and the 
required high temperature formation medium. Alternatively, iMolster et al.l (|l999[ ) have found evidence that 
the cooler, lower- density, disk-like material in some evolved stellar systems exhibits a high abundance of 
crystalline silicates at temperatures below the annealing temperatures, and have attributed their existence 
to an unidentified low-temperature crystallization process. Laboratory experiments have also indicated that 
partial crystallization can be achieved at room temperature through electron irradiatio n, although it is 
unclear whether this process is effective in astrophysical environments (jCarrez et al.ll2001l ). 



The third aspect which remains unclear is the iron content in the material. In addition to constraining 
the species and derived expected behavior of the quasar absorber dust material, it would be interesting to 
compare the iron content in silicate dust with the Fe/Si ratios inferred indirectly from gas phase deple- 
tions in other systems. (Although we note that some iron m a y be locked up in non-silicate dust, such as 
metallic or metallic oxide grains.) For instance, IWelty et al.l (|2001l ) have illustrated that along one sight 
line in the SMC, Fe exhibits severe depletion while Si can exhibit much milder depletions than would be 
expected from the depletions of other refractory species, in contrast with what is observed for the Milky 
Way. Galactic crystalline silicates are generally Mg-rich and Fe-poor, i.e. to be much closer to forsterite 
than to fayalite, although so me large (>1 /im diameter) cr ystalline silicates with F e are found in the Sola r 
System (jMolster et al.ll2002r ). However, many DLAs [e.g., iProchaska et al.l (1200711 : iKulkarni et al.l ( l2010l )] 
show at least some Fe-depletion, i.e. [Fe/Zn]~-0.62 for z<1.5 ([Meiring et al.ll2009l ). so their dust grains are 
unlikely to be completely Fe-poor. Longer wavelength IR spectral data are required to establish the relative 
Fe-enrichment in this source, as described in Sj5j 

Finally, there remain a number of alternatives for the formation mechanism of the crystalline silicates 
within the PKS 1830-211 absorption system. In general Galactic crystalline silicates are believed to form 
via one of three primary mechanisms: (i) condensation at or above the glass temperature, which is thought 
to be the primary mechanism around evolved stars; (ii) annealing (crystallization via heating), or formation 
following melting, of amorphous grains in a high temperature environment, which is the primary mechanism 
invoked in the accretion disks of young stars; or (iii) a by-product of planet-formation mechanisms, as invoked 
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for comets, such as through the coUisional cascade of asteroid-sized objects due t o gravitational effects, or 
through flash-heati ng by shocks stim ulated by protoplanet-disk tidal interactions ("Molst er fc KempeilEooi 
iMolster et al.,,2002 : Bouwman et al.i2003,: Harker fc Desch 2002.) . Although these mechanisms can success- 
fully produce sizable quantities of crystalline silicates, these crystallines are rapidly amorphized or otherwise 
destroyed by a range of mechanisms including sputtering, evapora.tion, cosmic-rays, grain-g r ain collisions. 



super nova shock waves and adsorption of Fe (jMolster fc Kempeij l2005t iMolster et al.l l2002t iTielens et al 



19981 ) . In the ISM, crystalline silicates may be am orphized within ~40-90 Myr, by, e.g., heavy ion cosmic 



rays (jKemper et al.ll2004L l2005HBringa et al.ll2007h 



In the case of ULIRGs, the re latively high ISM crystallinity may be explained by enhanced ongoing 
star-formation. ISpoon et al.l (j2006[ ) ascribe the high percentage of crystallines to recent merger-triggered 
star-formation, and postulate that the amorphization process lags the injection of dust which is driven 
by the intense star-formation. They posit that this crystalline-rich dust is originating in evolved massive 
st ars, such as re d supe rgiants, luminous blue variables, and Type II supernovae. Numerical simulations 



Kemper et al.l (j201ir) . however, suggest that even under the assumption of intense star- formation rates 



(1000 M0 yr~-^) and highly efficient dust production by supernovae, they are barely able to explain the 
6.5-13% crystallinity observed in ULIRGs, and with more realistic, lower, values for these rates, the addition 
of a secondary heating source, such as an AGN, is required to fully explain the observed crystallinity. 

The scenarios which appear the most probable to explain the observed crystalline silicates in the PKS 
1830-211 absorption system run a wide gamut, (i) The first scenario is an extreme burst of star-formation 
within the SW spiral arm, which is heavily obscured by both dust and molecular material. Alternatively, one 
of the PKS 1830-211 lines of sight could trace a path directly through an isolated, high stellar mass cluster 
in which strong outflows have transported the crystalline silicates into a slightly cooler region, wherein the 
crystalline material has not yet amorphized. While intense star-formation has been invoked to explain the 
crystallinity in ULIRGs, there is no evidence for intense star-formation within the spiral arms of the PKS 
1830-211 absorber galaxy. Studies have previously indicated that the z— 0.886 absorber h o st is a typical 
late-type spiral, based on its position on the TuUy-Fisher relationship [see, e.g.. IWinn et al.l (|2002l )]. While 
Rv is relatively high (iFalco et al.l 1 1999), there is a lso evidence that the dust-to-gas ratio in this system is 
consistent with the Galactic value ( Dai et al. 20061) . Furt hermore, the SW arm is noted to be signiflcantly 
rich in molecules and even displays unusual isotopic ratios ( Muller et al.ll20o"6 . 2011 ). (ii) The second scenario 
is intense, recent, AGN activity within the absorber galaxy which produced outflows to export crystalline 
silicates fro m the central torus, o r a strong flux of UV radiation to anneal the amorphous silicates, as has been 
invoked by iKemper et al.l (|201ll ) to help explain the ULIRG crystallinity. However, there is no published 
evidence of a strong AGN in the absorber galaxy, (iii) The crystalline material could have been abundantly 
produced in young stellar disks. However, in this scenario, the material should be hotter, and potentially 
detected in emission, rather than absorption. Furthermore, this mechanism would again require a massive 
(undetected) starburst to produce adequate quantities of crystalline silicates to dominate the profile, (iv) 
The remaining mechanism which could be invoked is that of shock-heated material, perhaps from a localized, 
extreme shock resulting from a transient event. Sh ock-heating mechanisms have been used to explain the 
crystalline silicates in comets ( Harker fc Deschll2002 ). although we note that very strong shocks could destroy 
not only the crystalline structure but also the dust grains themselves. The origin of the shocking mechanism 
in the absorbing galaxy remains unclear, however, and detailed simulations would be required to determine 
whether a large scale shock could produce sufficient temperatures to form the crystalline structures without 
destroying the dust grains. We regard scenario (i) as the most likely given the high molecular content and 
large Ry value of the absorber. 
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4.2. Amorphous Silicate Absorption Combined with Atomic/Molecular Features 

While we believe that crystalline silicates provide the simplest explanation for the observed substructure 
within the 10 /im absorption feature, we now explore a variety of scenarios to assess whether the observed 
features could be reproduced without crystalline silicates. The first alternative explanation which we explore 
is the possibility that the observed structure is the result of a combination of a broad absorption feature, 
produced by an amorphous silicate of olivine composition, and a series of other additional atomic and molec- 
ular absorption/emission features within the quasar absorber material. This possibility is worth exploring 
because (i) the 10 /im amorphous silicate feature is present in many G alactic and e xtragalactic sources; (ii) 
this a bsorption system is exceedingly rich in molecular material ( Wiklind fc Combesl[T996. ,1998 Muller et al 



20111 ): (iii) narrow absorption and emission features have been observed in ULIRGs dSpoon et al.ll2006l) : and 



(iv) previous detections of crystalline silicates in astrophysical sources, e.g. ULIRGs, have been much weaker 
than in our profile. 

The first test of this scenario which we have implemented is to manually mask potential absorption 
and emission features within the broader 10 /im absorption feature. We have then attempted to fit the 
remaining points with a series of four amorphous olivine profiles which were selected as providing the best 
range of olivine templates (see Appendix Table [TU] for template details). We have implemented these fits 
applying two different masking combinations: mask A which attributes the 9.3-10.05 /im and 10.35-10.85 /^m 
regions to unidentified, superposed emission features, and mask B, which includes both of these purported 
emission features, as well as an absorption feature at 10.85-11.3 /xm. (We have also considered two other 
masks which exclude the 8.65-9.25 /im spectral region, in which the enstatite, silica, or serpentine material, 
or ammonia, absorption, may be present, instead of the 9.3-10.05 /im region; but these two masks do not fit 
the data as well as mask A or mask B, and so we do not illustrate these fits.) We have adopted the same 
fitting procedure as utilized in These fits are illustrated in Figure [9] and parameterized in Table |6l We 
present the Xr values in this table in order to facilitate comparisons between the different templates and 
masking schemes, but we note that the bulk of the structure in the 10 /im region is masked, and as a result 
relatively few points in the masked profiles are being considered in the Xr- the less aggressive mask A, 
the amorphous olivine silicate fits remain too narrow/shallow relative to the PKS 1830-211 profile and peak 
at shorter wavelengths. 

As a second test we have employed the same masking of the PKS 1830-211 absorption spectrum, and we 
have additionally allowed the amorphous olivine templates to be arbitrarily shifted in wavelength-space to 
produce the best fit. We permit this wavelength-shifting because (i) as discussed in Appendix [Cl variations 
in silicate grain morphology can produce marked offsets in the peak of the optical depth profile; and (ii) 
recent work by R. Nikutta (private communication) suggests that sizable wavelength shifts in the 10 fj,m 
peak can be produced around quasars with a clumpy dust torus model which combines simultaneous silicate 
detections in absorption and emission, in conjunction with the rising quasar continuum. Although we are 
not looking at dust within a quasar torus, the work of Nikutta raises the possibility that a combination of 
absorption and emission by clumpy dust could contribute to the observed wavelength offsets. The results of 
this fitting are shown in Figure [10] and are tabulated in Table [6] We find that using mask B with shifting 
of +0.1 to -1-0.5 /im, we can achieve relatively low Xr values. However, as previously noted, this is based on 
substantial masking of the absorption feature, and there are few points which are being directly fitted within 
the deepest part of the optical depth profile; the bulk of the fitting leverage is being applied by the wings. 
Given the reduced number of sample points being fit, we do not consider the masked fits to be as reliable as 
those in our primary analysis. 



- 18 - 



Furthermore, the masking scenario is contingent upon ah of the identified emission and absorption 
features being physicahy exphcable by molecular and atomic transitions. In order to determine whether there 
are plausible atomic and molecular lines which could produce these features, we have considered prominent 
lines which are seen in other extragalactic and Galactic sources (see Figure [TT]). As is evident from this figure, 
none of the commonly observed features are broad enough to explain the significant substructure within our 
absorption feature. We note that H2 S(3) emission, which is commonly seen in ULIRGs, may be present in 
our z=0.886 absorber, producing the small hump in the larger substructure peak. This emission, however, 
cannot explain the entire substructure pe ak, as the line is no t only relativel y narrow in compa rison, based on 
the line widths in star- forming galaxies ([Smith et al.l 120071 ) and ULIRGs (jSpoon et al.l 120061 ) . but it is also 
not centered on the broad substructure peak seen in our data. We also do not see prominent H2 S(2) or H2 
S(l) emission at longer wavelengths, although these features could be lost in the noise. The emission feature 
from 10. 35-10.85 iim could be associated with the [S IV] 10.51 iJ,m fea ture, which is again seen in some 
ULIRGs (jSpoon et al.ll2006f ) and star- forming galaxies (jSmith et al.ll2007l ). but our feature appears broader, 
offset to longer wavelengths, and more prominent than the unresolved line seen in these other systems. The 
best candida te to explain the a bsorption feature from 10.85-11.3 /im remains crystalline silicates, based on 
the work of ( Spoon et al. 20061 ): however, we note that if we invoke a combination of amorphous silicates 
with molecular/atomic transitions to explain the bulk of the substructure in the profile, we require a much 
smaller contribution from crystalline silicates than claimed elsewhere in this paper. 

In general, these superposed emission features in conjunction with broad amorphous absorption appear 
inadequate to explain the detected substructure. Also, there is no strong evidence that all of these particular 
features can be physically produced and detected. Emission lines of II2 and [S IV] are detected in ULIRGs, 
but these objects are undergoing extreme star-formation, and can be detected in emission. By contrast, 
the PKS 1830-211 quasar absorber host is a regular spiral which is barely detected in HST images, with 
no noted pronounced star-formation. It seems unlikely that strong emission lines of these species would be 
produced in such an environment, and certainly not in the required line strengths/breadths to explain our 
substructure peaks. 

Lastly, we address the possibility that the substructure we attribute to crystalline silicates is, in fact, 
a by-product of PAH emission overlaid on a broad silicate absorption feature. Given the location of the 
expected PAH emission features in the rest frame of the z=0.886 absorber, illustrated in Appendix Figure [T5| 
we deem this unlikely. While the 10.68 /xm PAH emission feature in the rest frame of the absorber is consistent 
in location with the peak between the 10 /im and 11 /im absorption sub- features, this would still require 
a shift in the amorphous olivine template and an additional absorption mechanism to explain the '~9 /im 
feature. The sub-peak between 10 and 11 /im (in the z=0.886 absorber rest-frame) a hgns with the 5. 7 /im 
PAH feature in the quasar rest-frame, but this feature is quite weak in most quasars (jHao et al.ll2007l ). We 
deem it more likely that our substructure originates from crystalline silicates than from a wavelength-shifted 
amorphous olivine profile in combination with superposed PAH emission, although we concede that this 
sub-peak may be enhanced by some PAH contributions. 



4.3. Foreground Galaxy Absorption 



An alternative scenario to explain the structure in the PKS 1830-211 absorber spectrum is that we are 
detecting superposed amorphous silicate features from two absorption sources at different redshifts which 
serendipitously combine to produce the appearance of multi-peaked structure. This is p lausible in that th e 
NE quasar sightline is passing within 10 kpc of the center of a z=0.19 foreground galaxy ( Lovell et al. 1996I) . 
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and so our spectrum is also probing the material in the outskirts of this z=0.19 foreground galaxy, which 
may contain silicate dust. If we perform a bi-variate fit in which we consider an amorphous olivine template 
at the rest-frame of our z=0.886 quasar-absorption-host, combined with an amorphous olivine template at 
the rest-frame of the z=0.19 galaxy, as in Appendix IB. 21 we do not find this to be viable, however. At the 
wavelength at which we would expect 10 //m silicate absorption in the z=0.19 foreground galaxy, we instead 
detect an emission-type feature. We, therefore, conclude that the structure in the 10 fj,m region of the PKS 
1830-211 absorber spectrum is not the superposition of 10 fj,m silicate absorption from the absorber host 
galaxy and 18 /im silicate absorption from the z=0.19 foreground galaxy. We also note that while such a 
superposition might account for two of the substructure peaks, it would not account for three, and would 
require a third absorption system to explain the ~9 /xm absorption feature, or once again invoke additional 
mineral/molecular absorption in this region. Furthermore, while in principle two absorption features could 
combine to produce the observed 'W' shape formed by the 10-11 /ini substructures, individual amorphous 
olivine silicate profiles are too intrinsically broad to produce a 'W' of the required shape. 

We have also repeated the bi-variate fitting both allowing for a foreground galaxy at any arbitrary 
redshift, in combination with the z=0.886 absorber, and allowing for two systems at arbitrary redshifts out 
to z=l, with possibly no absorption from the z=0.886 system. These scenarios cover the possibility in which 
a dark object which is rich in silicate material was missed in previous studies, as well as the possibilit y that 
the z= 0.19 galaxy might in fact be at a higher redshift closer to the z=0.886 galaxy, as suggested by 



Lehar 



(|2000l ). or that the z=0.886 galaxy might in fact be multiple galaxies, as proposed bv lCourbin et al.l (|2002i ). 
However, we do not find any such "hidden" object responsible for producing superposed silicate absorption 
features in our spectrum. 

Finally, we have examined the possibility that we are detecting a combination of crystalline silicate 
absorption at z=0.19 with amorphous olivine absorption at z=0.886. While we do see structure near 6- 
7 /zm, which we have attributed in Appendix IA.3I to H2O and other molecular absorption, we note that 
these features occur in the z=0.19 rest-frame 10 /xm spectral region. We, thus, tested whether this structure 
might be attributable to crystalline silicates at z=0.19. However, when we fit the hortonolite spectrum to 
the z=0.19 rest-frame absorption spectrum, we find that the hortonolite features are too narrowly spaced to 
explain the ^6-7 /im features, and that the structure at longer wavelengths is not well-reproduced by 18 /im 
crystalline silicates in the rest-frame of the z=0.19 galaxy. Thus, we discount this scenario. 



4.4. Amorphous Silicates with Superposed Emission 



Lastly, we consider the possibility that the PKS 1830-211 spectrum could be produced by a broad 
amorphous olivine feature associated with the z=0.886 absorber in combination with superposed foreground 
or coeval emission. For instance, if the presumed z=0.19 galaxy is instead a massive galaxy located at 
z=0.886, as proposed by 
galaxies, as proposed by 



o r the z=0.886 host galaxy is instead comprised of two neighboring 



Courbin et all |2002 h. then we might expect a blend of silicate absorption with 



either some limited emission or absorption reflecting a differing chemical composition, from the second 
galaxy. This scenario has largely been excluded by the bi-variate fitting tests in ii3.3.4[ in which we did 
not find compelling evidence arguing for two distinct chemical compositions, although we cannot exclude a 
weak contribution from a second source. Alternatively, we could be detecting broad silicate absorption along 
the SW quasar line of sight probing the molecule-rich material, while simultaneously detecting emission 
from, e.g. crystalline silicates, along the other quasar line of sight. However, the broad amorphous silicate 
absorption would still need to be offset in wavelength to fit our spectrum, as would the crystalline silicate 
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emission template; and no single offset or explored crystalline silicate could simultaneously reproduce all of 
the observed peaks/substructure. Furthermore, as in the case of invoking atomic and molecular emission 
lines, it is not clear that such emission should be detectable from the z =0.886 host galaxy, or any companion 
galaxies. 

Alternatively, some of the foreground emission could be produced by a transient object within our 
own Galaxy or Solar System, perhaps connected to the scattered light feature discussed in Appendix IB. II 
The flux from a foreground object could be probing a longer rest-frame wavelength r egion than the z=0.886 
absorber, where crystalline silicate resonances would be more abundant. For instance, Bouwman et al. ( 2010l ) 



find abundant crystalline silicates around a young, low-mass M4 star with a proto-planetary disk, and as 
discussed in i 3B.2| there is an M4 dwarf within our slit. However, we do not believe that the flux from this 
particular star contributes significantly to our PKS 1830-211 absorption spectrum, as detailed in §B.2| nor 
is it believed to be a particularly young stellar system, and so any crystalline silicates produced during its 
youth are expected to have dissipated. Furthermore, the likelihood of a transient source providing significant 
crystalline structure seems unlikely, and as discussed in Appendix lB.il we do not expect contamination from 
the scattered light feature to be strong enough to produce structure at the level we observe in the PKS 
1830-211 spectrum. 



SUMMARY 



In summary, based on an infrared spectrum of PKS 1830-211 obtained with the Spitzer IRS, we have 
detected the presence of a 10 /zm silicate absorption feature in the z=0.886 absorber at a > 10— cr significance 
level, as well as the presence of an 18 fim silicate absorption feature. This reinforces previous findings from 
Kulkarni et al] |2007bl 120111 ) that silicate material is a component of dust in quasar absorption systems. We 
have found that the 10 /zm silicate feature is best fit by crystalline olivine, possibly in combination with a 
second material such as enstatite, silica, or serpentine, or with superposed ammonia absorption. Given the 
rarity of such highly crystalline silicate dust, we have explored a variety of other scenarios invoking atomic and 
molecular absorption/emission lines, and a superposition of silicate dust in two quasar absorption systems; 
but none of these alternative scenarios satisfactorily reproduce the observed structure in the PKS 1830-211 
spectrum. 

In order to confirm the presence of the postulated crystalline silicates, and to eliminate some of these 
suggested formation scenarios, supplemental data are required. In particular, moderately high resolution 
data of the PKS 1830-211 absorption system at >20 /im, in the z=0.886 absorber rest frame, are needed to 
provide evidence for the presence of longer wavelengt h crystalline resonances and to constrain the relative 
amount of Fe present in the system. As illustrated in jjager et al. I Il998) . at >20 iiui the crystalline olivine 
profiles vary significantly with the ratio of Mg/Fe, with peaks shifting to longer wavelengths , by a s much 
as several /im, as more Fc is introduced to the mineral. As summarized in lMolster fc Kempen (j2005[ ). these 
differences have been used to provide evidence of the hig h Mg/Fe ratio in, e.g., evolved stars, using the 69 /im 
band. Furthermore, as discussed in iKoike et all |2006l) . variations in the temperature of the material will 



produce offsets in the peak wavele ngths of the spectral features; these offsets become increasingly sizable at 
longer wavelengths. For example, iBowev et al.l (|2002al ) have exploited these offsets to place constraints on 
the temperatures of a range of stellar objects using the 69 fj,m. forsterite band. Astronomical observations 
spanning a broader spectral range could also provide evidence to support/rule out other atomic, PAH, and 
molecular transitions in both the quasar and the absorber. Higher resolution spectra covering the absorber 
rest-frame 10 /im and 18 /im spectral features for the PKS 1830-211 z=0.886 absorber could also be of 
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great use for ascertaining better estimates of the structure features and line widths, which in turn could 
help discriminate between the emission/absorption scenario and the crystalline silicate scenario, as well as 
helping to determine whether we we are seeing broad emission line features, or emission features which 
are a blend of multiple narrower emission lines. These data would also help rule out any contributions to 
the spectra from the anomalous scattered light feature (Appendix IB.1[) . as it is believed to result from a 
transient source, and so should be absent in follow-up observations. Finally, spectroscopic or photometric 
observational constraints on the amount of star-formation within the obscured SW spiral arm of the z=0.886 
absorbing galaxy could help to further constrain the formation mechanisms of the silicate dust. 
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A. PKS 1830-211 Quasar-Continuum Normalization 

A.l. Spectral Joining Methodology 

As addressed in §2.2| we have joined the four individual nod-combined SL and LL spectra in order 
to form a single spectrum for PKS 1830-211, following a manual trimming of the data to eliminate lower 
quality data points, i.e., those points which appear visually discrepant or noisy. As discussed further in 
§A.4[ the exclusion of these data points does not significantly impact our derived fits and conclusions. The 
majority of the removed points constitute isolated data points where a bad pixel was not removed through our 
conservative spectral cleaning described in ^2.l\ The remaining points occur in regions where the two nods 
showed discrepant spectral features, indicative of bad-pixel contamination in one position, or occur beyond 
38 /im where the noise notably increases. Lastly, we have eliminated points at the edges of each spectral 
order which do not connect smoothly with the adjoining orders. As in ^2.1[ we have been as conservative as 
possible when making these cuts, so as not to eliminate spectral features which may be potentially physical 
in origin. We illustrate the spectrum which would result when excluding the bulk of this clipping in the left, 
bottom two panels of Figure [12] (clip v3 and v4). As is evident in this figure, there are several additional 
data points with large error bars, which are absent in our fiducial clipping (vl), but the overall shape of the 
spectrum is unaltered. 

Following this manual clipping of the data, we have normalized and stitched together the data in the 
SL2, SLl, LL2, and LLl spectral orders, to form a single SL spectrum and a single LL spectrum for the PKS 
1830-211 system. Within the SL and LL orders we have normalized the first and second order data to the 
bonus (third) order. The individual data points from the bonus orders have been excluded from our final 
spectrum, as they provide no additional information since they align with the data in the retained spectral 
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orders and exhibit similar features within the overlapping spectral regions; furthermore, their inclusion 
would lead to uneven sampling over the spectral regions spanned by the bonus order relative to the full 
spectrum. We note that the joining of the LLl and LL2 spectral orders occurs within the 10 fim feature 
of our z=:0.886 absorber, which is the focus of our analysis, and so could potentially impact the shape of 
this feature. However, as the orders align well, have overlapping data points at the edges of the orders, and 
match the bonus order spectrum, we do not believe that the observed substructure illustrated in Figure [l] 
and attributed in S|4]to physical mechanisms, is affected by this joining. 

We have generated the final PKS 1830-211 spectrum by scaling down the LL spectrum and joining it with 
the SL spectrum. An examination of the well-separated longer (radio) and short (X-ray) wavelength data 
for PKS 1830-211 in the literature does not convincingly argue for either scaling up the SL data or scaling 
down the LL data. Since the SL data generally have smaller associated formal uncertainties, we have given 
that order preference in setting the absolute-flux-normalization. This has resulted in final normalization- 
scaling factors of 0.97, 0.99, 0.91, and 0.87 relative to the original extracted data for the SL2, SLl, LL2, and 
LLl data respectively (including the shifts of the first and second order data to match the corresponding 
bonus order data), with a slight overlap in SL and LL data points at th e location of the jo i ning. These 
scalin g factors are consistent with scalings employed in other analyses [e.g.. Spoon et al. ( 20061) : Smith et al 



(|2007i )]. As our analysis is based on a quasar-continuum-normalized spectrum, any offsets in the absolute 
flux normalization will not impact our conclusions. We have also considered a slightly modified version of 
the joining in which we employed additional clipping at the long-wavelength edge of the SL spectrum, and 
then shifted the LL spectrum down by a larger factor. This alternative joining is illustrated in the clip v2 
(and clip v4 for the less aggressively clipped data) panels of Figure [121 ^nd does not significantly impact our 
conclusions. 

We have, lastly, experimented with binning the data by a factor of three, using an unweighted average. 
While, as illustrated in the right panel of Figure [12] this can reduce the uncertainty associated with each 
data point, and produce a cleaner looking spectrum, it can also obscure the finer structure in the absorption 
feature. Given that we find indications of significant substructure in our analysis, we have opted against 
binning. As discussed in ^^i^ however, had we chosen binning it would not have altered our conclusions 
that crystalline olivine provides the best match to the 10 /xm data feature. 



A. 2. Variations in Spectral Joining and Continuum Normalization 

In order to fit the silicate template profiles in i j3.1l to the PKS 1830-211 absorber silicate features (in the 
rest frame of the z=0.886 absorber), we have normalized the observed spectrum by a fit to the underlying 
continuum. We have first examined the data at both longer (radio) and shorter (X-ray) wavelengths to as- 
certain whether a single quasar continuum could be fit through an extended frequency range. However, given 
the discrepancies in the literature, the known temporal variability across a wide spectral range of the quasar 
and u ncertainty in the shape of the quasar spectral energy distribution in the IR region (i Ponnarumma et al 



201 II ). as well as the large wavelength separation between our data and these supplementary data, we decided 



to select a continuum which fits in detail solely our own IR spectral data. 



We then explored a range of more than 20 fits to the quasar continuum. In the analysis of lKulkarni et al 



(|2007bL 120111 ). the underlying quasar continuum, in the quasar absorption systems, was well- represented by 
a polynomial or power-law. Thus, in this analysis we considered not only a linear polynomial and a power- 
law normalization, but also 3rd-5th order Chebyshev polynomials, 3rd order Legendre polynomials, and 
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single component cubic spline fits. The power-law and linear polynomial fits were determined using a least- 
squares algorithm, while the higher order polynomial fits were obtained using the IRAF task continuum with 
unweighted fitting in linear space, and no averaging of the sample points. A subset of the most representative 
of these fits is tabulated in Table [7] and illustrated in Figure [T31 Unless specifically noted, each of these fits 
was obtained using masking to exclude the four prominent absorption features in the spectrum: 5.5-6.2 /im, 
6.5-7.1 /zm, 8.5-12.2 /im, and 14.8-18.3 fj,m in the z=0.886 absorber rest-frame. We also forced the IRAF- 
fitted curves to pass through points longward of 18.5 /im in order to constrain the curvature of the continuum 
normalization; the implications of this constraint are addressed below. Based on the relative contributions 
of the continua to the observed and laboratory fits over solely the 10 /zm feature and over the combined 
10 and 18 features, as discussed in and a visual inspection of the resultant normalizations, we 

have determined that the 3rd order Chebyshev polynomial provides the optimal fit. As addressed in HAAl 
however, adopting an alternative normalization would not have altered our conclusions in this paper. 

Lastly, we have experimented with forcing the 3rd order Chebyshev polynomial to pass through/near 
certain critical parts of the spectrum. Since the continuum fitting is unweighted, this prevents anomalous, 
high uncertainty points, like those at the longest wavelengths, from skewing the fit. In Figure[T3]we illustrate 
five alternative assignments of these regions. The fiducial fitting solely forces the curve to pass through 
points at the longest wavelengths (A >18.5 /im), termed 'A' in the figure. We have also examined forcing the 
spectrum through points between the 10 and 18 features (termed 'B' in the figure), points just shortward 
of the 10 /im absorption feature (termed 'C'), points just shortward of the 6 H2O absorption feature 
(termed 'D'), and points in the substructure peak near 9.6 /im (termed 'E'). These variations are illustrated 
by the shaded grey regions in Figure 1141 9'iid listed in Table [7l with the resultant continuum variations 
shown by the red lines in the left side of Figure 1141 These variations can impact the relative depth of our 
absorption feature, as addressed in ^AAl and illustrated in the right side of Figure [TU but they do not alter 
our conclusion that the absorption feature in the PKS 1830-211 spectrum is best-represented by crystalline 
olivine. 



A. 3. Potential Residual 'Emission/ Absorption Bumps in Normalized Spectrum 



An examination of our final normalized spectrum, shown in the right side of Figure (TJ illustrates that in 
addition to the prominent 10 nm and 18 fim. absorption features, which are the focus of this analysis, there 
are several additional emission and absorption features, which we briefly address here. First, there are two 
absorption features located near 6 /im and 7 /im, respectively. The first of these we attribute to water ice 
(H2O). Galactic observations of, e.g., ice-rich embedded young stellar objects (YSOs), have also found this 
feature, and a direct comparison of these Galactic profiles with our absorption feature finds good agreement 
in both the breadth and shape of the H2O feature ( Gibb et al. 200d Chiar et al. 2011 ). The slightly longer 
wavelength feature is al so seen p rominently in some Galactic environments, alongside the H2O feature, e.g., 
in the W33A spectrum (jGibb et al. 2000). Although its origin is not completely understood in the literature, 
it is likely attributable to other molecular transitions associated with the ices. Given that the PKS 1830-211 
z=0.886 absorber is well-documented to be rich in molecules, we find that these associations are plausible. 

We also note the presence of a visible emission- type feature just shortward of the 10 /im absorption fea- 
ture, near 8 /im, which we contend may be evidence of PAH emission. As illustrated in Figure[T5l the location 
of the feature is consistent with the 7.7 /xm complex and 8. 3-8.8 /xm PAH emiss ion fe atures in the rest-frame 



of the z=0.886 absorber. Following the classifications of iPeeters et al.l (j2002l) and Ivan Diedenhoven et al 



(j2004l ) , the shape of the feature is similar to a blend of perhaps Type A and Type C PAH emission, perhaps 
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originating in a range of stellar sources. We note that if this feature does originate in PAH emission, it 
is unusual that we do not see evidence of the usually prominent 11.3 /im complex emission, although the 
relat ive strengths of the PAH features can vary, for example, due to variations in the ionization levels [see. 



Tielenj |2008[) ; ISmith et al.l (|2007l )] . These PAH ratio variations are illustrated in Figure 1161 comparing 



e.g., 

the PAH features for two star-forming galaxies, NGC 5866 and NGC 3198, rescaled and overlaid on top of 
our PKS 1830-211 spectrum. Furthermore, it is possible that our absorption feature could be hiding some 
weak PAH emission, i.e., that the rightmost (^11 fim) substructure peak may indeed be deeper or broader 
than our spectrum indicates. This would not alter our assessment that the primary constituent of the sil- 
icate dust is crystalline olivine, but might alter our ascribed species of crystalline olivine or the requisite 
se condary materials (jSI ). Finally, we note that in the Q0852-I-3435 quasar absorption spectrum examined 



bv lKulkarni et al.l ()2011l ). there is also evidence of not only the two absorption features near 6 and 7 /xm, but 



also the emission features shortward and longward of the 10 fim absorption; thus, these emission features 
may be physically connected to an enriched environment. 

Finally, we address the 'emission' feature longward of our 10 /im absorption. The physical origin of 
this feature is less clear. It may originate from one or a combination of mechanisms including weaker PAH 
emission features in the rest frame of the absorber, emission from the 7.7 PAH complex emission in the 
frame of the quasar host galaxy, H2 S(2) emission in the frame of the z=0.886 absorber (which would be 
consistent with evidence for the H2 S(3) line at shorter wavelengths), or it may be related to the scattered light 
feature addressed in Appendix lB.il which provides the greatest contamination in this region, as is evidenced 
by the larger associated error bars. Given this range of plausible mechanisms, and the added complication 
of the contaminating scattered light, we do not attribute this feature to an artifact of the quasar continuum 
normalization. We note, however, that treating this emission as an artifact, and eliminating it through an 
altered continuum shape, is encompassed by our tests explored in Figure [HI and Table [71 and we emphasize 
that these normalizations would not impact our qualitative conclusions in this paper. 



A. 4. Impact of Data Treatment on Derived Peak Optical Depth 

Lastly, in order to ensure that our data processing choices regarding clipping, binning, and joining 
(Figure [T^. the quasar continuum polynomial (Figure [T5)) . and the continuum shape (Figure [H)) have not 
impacted our conclusions, we have repeated our template fitting described in using the alternative 
normalizations in Table [Tj We find that, as detailed below, while the continuum variations can have some 
impact on the optical depth, even the most extreme plausible normalizations do not alter our conclusions 
that the observed structure is best represented by crystalline olivine. 

We find that neither the possible clipping, joining, and binning options, nor the variations in the shape 
of the adopted polynomial normalization, significantly impact our peak optical depth measurements. All 
explored clipping, joining, and binning options increased the derived optical depth by values within one- 
sigma, and over the 10 fim region hortonolite remains the best fit. The linear and power-law normalizations 
produce optical depth measurements which are consistent over the 10 /im range, but we discount the fitting 
over the full fitting range because of their visibly poor fits at the longest and shortest wavelengths. We 
also discount the 5th order Chebyshev, as we believe that it has too much curvature and is over-fitting the 
continuum, and thus removing some of the residual substructure. The remaining continuum normalization 
fits are all consistent within one-sigma of the fiducial fit, and also consistently find that hortonolite produces 
the best fit. 
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Finally, we examine the impact of the different shape constraints on the third order polynomial. We 
find that the more aggressive forcing, which eliminates the emission peaks preceding and following the 10 /im 
absorption feature, can raise the optical depth to at most 0.42±0.05 using the hortonolite template, or to 
a maximum of 0.37 using the crystalline olivine template composition which provides the best fit. The 
Ar=0.10 is nearly consistent within the one-sigma uncertainties, relative to the fiducial fit. We also note 
that imposing these constraints we consistently find that crystalline olivine provides the best fit, although 
the precise species of olivine varies. We prefer our fiducial normalization to these alternatives, because it is 
able to better fit the lowest wavelength data in our spectrum, which is under-fit in these more aggressively- 
forced polynomials. If we employ less aggressive forcing, as explored in test 'fl' (see Table [TH)) . the 10 ^m 
optical depth (rio) decreases slightly to 0.26, while the fit including longer wavelengths (t/„;/) increases to 
0.34; it is at this long wavelength end that the continuum normalization differences are the most significant. 
We believe that our adopted polynomial, however, provides a more reasonable fit to the overall continuum, 
and is less influenced by noisy points at the longest wavelengths. Lastly, using a polynomial shape which 
passes through the ^9.6 /im substructure peak, we find that the optical depth is decreased to rio=0.20 and 
''"/ti(;=0.24, both of which are consistent within one-sigma of our fiducial fits, and again, argue for crystalline 
olivines as the origin of the substructure. 



B. Minor Caveats 



We note the presence of two distinct spectral contaminants, one physical and one spurious, which impact 
our data but which we believe do not significantly alter our conclusions. 



B.l. Scattered Light Feature 

First, we have identified an anomalous emission feature in our 2-D LL spectra which we associate 
with scattered light. The relatively broad feature is adjacent to our quasar spectrum at approximately 
(ct, i5) j2ooo=(18:33:36, -21:03:41.3) in each of our LLl, LL2, and LL3 spectra. The centroid of this feature is 
approximately 8 pixels away from the quasar spectrum at the wavelength relevant to our 10 fiia absorption 
feature. This prominent anomalous feature exhibits a spatially narrower distribution at low LL wavelengths, 
which broadens progressively as the wavelength increases; at the longest LLl wavelengths it appears to be 
bifurcated with the light generally concentrated in two spatially distinct prongs, and a diffuse contribution in- 
between. Thus, while at the shortest LL wavelengths we do not believe this feature substantially contributes 
to our spectrum, at the longest wavelengths it visibly crosses through our narrow PKS 1830-211 spectrum. 

We have determined that the origin of this feature is most likely scattered light, from a bright off-field 
object. A visual inspection of both ground-based and HST visible and near-IR data reveals no significant 
objects at the spatial location of the feature. Furthermore, we have utilized SPICE to extract the spectrum 
of this feature, an d compared it with a catalog of known stellar spectra observed with the Spitzer IRS 



(jArdila et al.ll2010[ ). and with quasar spectra; the feature does not match any of these spectral shapes. The 
Spitzer Science Center help desk has confirmed that this feature is in all probability a very rare scattered 
light feature, of unidentified origin, although the most likely candidates are either an IR bright source near 
the galactic plane, or a transitory object from our own solar system. Neither we, nor the help desk, could 
find any bright optical or IR sources in the immediate vicinity of our target. We have also ruled out detector 
features, such as latents or bad dark corrections, as the origin of this feature. 
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As this scattered light feature is spatiahy broad, and visibly crosses through the quasar spectrum at the 
longest wavelengths, it may be the origin of the observed discrepancies between the two LLl nod positions 
in the 22-26 /im spectral region. This location corresponds to the 11.7-13.8 fj,m region in the rest frame 
of the z=:0.886 absorber. As discussed in W2.21 a discrepancy on the order of 15-20%, in which one nod 
position exhibits a broad emission-type feature while the other exhibits a broad absorption-type feature, is 
present over this spectral range. In order to identify the origin of this discrepancy, which is larger than 
those observed in the SL and LL2 spectral orders, we have first examined the position of the slit-centering to 
determine whether perhaps an offset in centering is causing a foreground/background object to fall within 
the slit for one of the nod positions, but not for the other. However, the header information indicates that 
there are variations of at most (Aa,A(5)=(0."018, 0."009) and (0."020, 0."014) within the two nod positions, 
and no significant offsets (within one pixel) in-between the nod positions; IRS exhibits pointing uncertainties 
on the order of 1" (which is smaller than the pixel scale). There is no evidence to suggest that we have a 
slight pointing mismatch, which could result in the inclusion of extra flux from a foreground/background 
object in only some of the exposures, nods, or orders. 

We have also explored whether the difference between the two nod positions could stem from the 
background subtraction. As we were unable to subtract the background by differencing the two nod positions 
because of the anomalous scattered light feature, we have instead constructed a background generated from 
the combined off-source LL2 images. It is plausible that a background gradient, independent from the 
scattered light feature, might be contaminating one nod background or region of the chip but not the other. 
In order to check this, we have utilized SPICE to extract a background spectrum from multiple locations 
across our LL2 background image, as well as across a region of each of the LLl spectra not covered by 
the quasar or the scattered light feature. We do not detect any substantial variations between the intrinsic 
backgrounds in the two nod positions, nor do we detect any significant spatial variations in the subtracted 
background, indicative of either a gradient or of a region of slightly bad pixels underlying one of the nod 
positions. We have also confirmed that the feature is not a result of any treatment of the background by 
the extraction algorithm, by replacing the pixels within the scattered light feature of the 2D spectrum with 
representative background pixels, and repeating the extractions; the extracted structural features remain 
unchanged. An additional test using the IRAF task apall confirms that treatment of the background in our 
extraction algorithm is not the origin of the derived structure. 

Additionally, we have explored whether the difference in the 22-26 /zm spectral region could stem from 
a localized event along one of the spectra, such as an isolated cosmic ray burst or transitory event. In order 
to test this, we have extracted the combined 20-exposure spectrum, for each nod, as well as four subset 
spectra with 5 exposures each. We find that all of the spectral features in each nod position, including the 
discrepancy, are virtually identical between these five extractions. Therefore, it is unlikely that an isolated 
transitory event produced the discrepancy. 

Thus, the most likely explanation for the discrepancy in the 22-26 /Ltm region is contamination from 
the scattered light feature, which peaks in intensity over this 4 /xm region, and more strongly affects one 
of the nod positions than the other. An inspection of the 2-dimensional spectrum shows that the scattered 
light feature is both the most narrow and has the highest surface brightness at the short wavelength end 
of LLl, where it is well-separated (~8 pixels) from the quasar spectrum. As the wavelength increases, the 
breadth of this feature also increases. However, the potential impact of this increased contamination at 
the longest wavelengths is mitigated by the decrease in surface brightness of the feature with increasing 
wavelength. In order to determine whether a narrower extraction window might minimize this contribution, 
we have compared a SPICE extraction obtained using the default point-source extraction template, which 
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is of variable spatial width becoming broader at longer wavelengths, with two manual template extractions, 
one of which is 3 pixels wide, and one of which is 5 pixels wide. Modifying the extraction template impacts 
the absolute calibration of the extracted spectrum, but as we are only utilizing the manual extractions to 
compare features within the spectrum, this is not an issue. While we do find that the continuum level is 
slightly different between the extractions, as expected, we find no significant differences in the shape of the 
PKS 1830-211 spectrum or in any of the spectral features, along either nod position. Thus, employing a 
narrower extraction template would not change our analysis. We also have no method by which to separate 
any possible contamination occurring directly over the location of the quasar spectrum from the quasar 
spectrum itself, as this scattered light feature does not exhibit any standard spectral type. We, thus, rely 
on the combination of the two nod positions to effectively cancel out the contamination pattern within this 
22-26 fim region, with resultant larger uncertainties associated with the data points. 



B.2. Contamination from Foreground Objects 



The second contaminant of our quasar absorption system spectrum comes from several foreground ob- 
jects which are also present within the slit. However, none of these objects seem to be bright enough to 
affect our analysis, as we describe below. As discussed briefiy in i )2.11 in addition to the two images of 
the lensed quasar and the z=0.886 quasar absorber host galaxy within the slit, a foreground star and a 
lower-redshift foreground galaxy are also encompassed by the IRS slit in all spectral orders. The brighter 
of the two quasar sight lines, termed the NE or A component in the literature, is located at 18h33m39.931s 
and -21°03'39."75 (J2000), while the dimmer, and more dust- obscured quasar sightline to the SW is lo- 
cated at (Ap, A^ ) = (-0."64, -0."73) from the NE component |Subrahmanvan et al.l[l99ol : Ijin et al.lllQQQl: 
Winn et al.l 120021 ). The host galaxy for the absorber, which has been ide ntified as an Sb/S c spiral, is sep- 
arated by (Aa, A(5) = (-0."33, -0."49) from the NE quasar component (jwinn et al.ll2002h . The absorber 
subtends a physical scale of 7.6 kpc/", assuming Ho=73 km s^^ Mpc^^ and a standard ACDM cosmology. 
Our requested central slit position of 18h33m39.89s, -21d03m40.4s places it close to the center of this quasar 
absorption system. 

The foreground star which is also present within the slit, termed SI in the literature, dominates the 
combined light distribution at optical wavelengths, but has a minimal impact at the mid-IR wavelengths 
under consideration in this analysis. This star, which is located at (A a, A^) = (0."09 , 0."53 ) rel ative to the 



NE qu asar component, has been identified as a likely M4 dwarf star bv lCourbin et al.l (j2002l ) and lDiorgovski 
(|l992l ). Star SI dominates the total light at visible wavelengths, contributing between 54-74% of the total 
light at F555W (V) and 68-83% of the total light a t F814 W ( I), based on th e mag nitudes for the star and 
other objects within the slit provided bv lWinn et al.l (|2002l ) and lCourbin et al.l (|2002| ): differences in the total 
contribution by SI reflect variations in separating the light contributed by different objects by these authors. 
In the near-infrared, however, SI contributes a much smaller perce ntage of the total light: 32-44% at F160W 
(H) and 17-23% at F205W (K) (jLehar 2000: ICourbin et aLlliooj ). Extrapolating these rapidly decreasing 
relative contributions to the spectral region covered by the IRS, we would expect ^1% of the total fiux to 
be contributed by SI at the shortest (5/xm) wavelengths covered, and <0.01% at the longest wavelengths 
covered. These findings of a negligible contribution to the total flux are reconfirmed by rescaling the IRS 
spectra for three brighte r early M dwarf sta rs (IID180617, GJ687, and GJ849) provided in the compilation 
of IRS stellar spectra bv lArdila et al.l (|2010l ): we find wO.5% contribution to the total flux near 5 /im, and 
0.001% contribution to the total fiux near 35 /im. In the region where the z=0.886 quasar absorber 10 /im 
silicate feature is located, we find only an 0.01% contribution to the spectrum by these rescaled, similar stars 
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from the Spectral Atlas. The small contribution of the M4 dwarf to our combined spectrum is a result of 
the rapidly decreasing strength of the stellar emission at mid-IR wavelengths, with the increasing strength 
of the underlying quasar continuum flux density. Furthermore, we have verified through an inspection of the 
HD180617 stellar spectrum that we do not expect any significant structural features to be produced in our 
spectrum from an M4 dwarf star. We, thus, exclude the possibility that the spectral features we see in the 
PKS 1830-211 spectrum are produced by this foreground star. 

The secon d potential foreground contaminant to our q uasar spectrum is the z=0.19 iLovell et al.lllQQd 
Wiklind k Combes 1995 ICourbin et"al]|2002t Iwinn et alll2002 ) galaxy, termed G2 in the literature, which 
also lies within the slit, and is separated by (Aa, AJ) = (0."24, —2. "49) from the NE quasar component. 
We do not believe that this galaxy will produce a significant contribution for three reasons. First, it is 
substantially (Ami?205VK — 3.35 mag) fainter than the NE quasar sightline in the near-infrared, and is less 
reddened than the quasar (Lehar 20OT). Hence, any mid-IR emission produced by the galaxy should be 
negligible. Second, any spectral features contributed by the galaxy will be well-separated spectrally from 
those of the 0.886 quasar absorber host galaxy. Third, the central regions of the foreground galaxy are not 
being directly illuminated by either of the quasar sight lines, and so material from the foreground galaxy 
should not be producing significant absorption in the quasar spectrum. However, the closer, NE quasar line 
of sight passes within 10 kpc of the center of this foreground galaxy (assuming a standard ACDM cosmology 
with Ho=73 km s~^ Mpc~^) and so may be probing the outer regions of this z=0.19 foreground galaxy. Thus, 
in order to confirm that none of the prominent features which we observe coincide with silicate absorption 
in a z=0.19 galaxy, we have shifted the intermediate amorphous olivine profile (AmOliv) into the rest frame 
of the z=0.19 galaxy and repeated our fitting procedure. As illustrated in Figure [ITl we find that while 
the purported 10 /im feature could be explained by 18 /xm silicate absorption in the z=0.19 foreground 
galaxy, there is no corresponding 10 /im absorption feature for this foreground galaxy. At the wavelength 
at which we would expect to detect 10 /im silicate absorption in the z=0.19 foreground galaxy there is a 
weak emission-type feature detected instead. We, therefore, conclude that the structure we see in the PKS 
1830-211 spectrum does not originate from a superposition of 10 /im absorption from the quasar absorber 
and 18 /im absorption from the z=0.19 foreground galaxy. 

The deblending of the objects in this field, based on their HST WFPC2 and NICMOS i mages, remains 



uncert a in, and it is possible that t here are furth er additional objects contained within the slit. ICourbin et al 



|2002h . iMevlan et all |2005l ) and iLehail |2000l ) lave identified a second for eground star terrn ed star P, lo- 
cated at (Aa, /S.5) = (— 0."3, — 0."5) relative to the NE quasar component; I Winn et al.l (|2002l ) have instead 
argued that this object is the bulge of the z=0.886 galaxy, and the current photometry is unable to distin- 
guish between the two scenarios. Given that this object is 0.02-0.05 times fainter than star SI at F205W, 
its emission would not have a sig nifican t impact on the derived profile should it indeed be a foreground 



star. Furthermore , ICourbin et al.l (|2002l ) argue that the object considered to be a single z=0.886 galaxy by 
Winn et al.l (|2002r ). may in fact be comprised of two galaxies near z=0.89, one of which exhibits significant 
spiral features. The impact of two companion galaxies is addressed briefly in 21 in which we have searched 
for simultaneous absorption, or simultaneous emission and absorption, at arbitrary redshifts in order to ex- 
plain the observed spectral features. However, this explanation does not adequately account for the spectral 
structure which we observe, which indicates that if two galaxies are contributing to the observed silicate 
absorption profile, they either have a similar chemical composition and redshift, or one of the objects is 
dominating the absorption profile. 
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C. Expanded Selection of Template Fits 

In this appendix, we tiave considered an expanded set of observationally-derived and laboratory optical 
depth templates. While this expanded analysis does not alter our fundamental conclusion that crystalline 
olivine provides the most viable explanation for the structure in the 10 /im silicate absorption feature, it 
allows us to explore how variations in environment, temperature, and silicate grain shape, size, and metallicity 
can impact our fits. We have considered this large selection of template profiles to encompass the broadest 
range of both the physical environments in which silicates are detected, and of the observed silicate chemical 
compositions. The degree of crystallinity in these silicates varies widely amongst different astrophysical 
sources, as discussed in detail in ^ ranging from nearly 90% for some comets to nearly 0% for the ISM. 
Thus, our broad range of template profiles implicitly allows us to probe a range of crystalline:amorphous 
silicate ratios, although we do not expect that the formation mechanisms producing the silicate absorption 
towards PKS 1830-211 are necessarily similar to those in, e.g., a comet. Reflecting the range of chemical 
compositions and grain properties producing the putative 10 /xm silicate feature, the observed profiles span a 
range of values for the peak feature wavelength from 9.37 /im for a YSO embedded in the Taurus molecular 
cloud to 10.20 fim for a ULIRG. The laboratory silicate optical depth profiles show an even broader range 
of peaks. We note that in our PKS 1830-211 spectrum, as discussed in the following sections, the peak 
wavelength of the absorption feature is longward of 9.7 /im. Laboratory experiments have found that 
variations in grain morphology, including size, shape, porosity, agglomeration, Si04 polymerization, and 
metallicity, can subst a ntially al t er bot h the shape of the optical depth profile, and its peak wavelength [e.g., 
Henning et al] |2005l) : iHennh^ |2010l) ]. 



C.l. Expanded Subset of Observational Templates 



We first examine an expanded set of observationally-based templates for astrophysical sources ranging 
from comets to star-forming galaxies. Although these profiles are inherently subject to more sources of 
uncertainty than the laboratory profiles, which are the product of controlled experiments, they can provide 
valuable insight into the astrophysical environment which might be lacking should we solely use laboratory 
templates. We list our expanded set of observational templates, separated by environment, in Table [5J 
with the derived fits listed in Table [9l and illustrated in Figure 1181 It is apparent both visually, and from 
examining the Xr of the fits, that none of these objects matches the profile in our system as well as crystalline 
olivine. This is significant because it suggests that our system is dissimilar from the Galactic ISM and from 
the ISM in several star-forming galaxies. 

The profiles which best match our system are the two comets and the ULIRG 060301-7934 (when 
excluding the absent 11.3 nm PAH feature). What unites these systems is that they are all crystalline-rich. 
While we do not have a precise estimate of the crystallinity for the Bradfield and Levy comets, comets, in 
general, are found to be up to 90% crystalline. The ULIRG 060301 -7934 has a 13% cry stallinity and provides 
a better fit than either ULIRG 18443-^7433 or ULIRG 00397-1312 |Spoon et al.ll2006h . both of which exhibit 
slig htly lower levels of crystallinity. These three ULIRGS were selected as the best matches from those in 



the Spoon et al. ( 20061 ) study; the ULIRGs with 6.5-9% crystallinity generally provided even poorer fits. 



The Galactic and stellar environments generally provide relatively poor fits. The best fit is from the 
AGB star, and since AGB outfiows have been found to be crystalline rich (10-15%), this is consistent with 
our need for silicate crystallinity. The diffuse ISM sources, wh ich were adequate for reasonably fitting the 
silicate absorption features in the iKulkarni et al.l (|2007bL 120111 ) quasar absorbers, are a poor fit here, being 
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shallower and offset to lower wavelengths relative to our feature. The circumstellar environments, dense 
molecular clouds, and stars embedded within molecular clouds all suffer from this s ame general failin g. 
Some of the molecular clouds, such as p Oph, characterized as a disturbed environment ( Bowev et al.ll2003i) . 
and the visually similar (not depicted) Serpens profile, produce significantly inferior fits when compared 
with, e.g., profiles for sources probing the the Taurus molecular cloud. Interestingly, profiles for sources 
embedded within the dense molecular clouds tend to produce slightly better fits than those illuminating 
the dark molecular clouds, which may stem from the fact th at young sources su r round ed by circumstellar 



material can be rich in crystalline silicates [e.g., see review in iMolster fc Kempeil (|2005[ l] 



C.2. Variations in Laboratory Amorphous Olivine Templates 

We also explore, in the context of the amorphous olivine silicate profiles, the range in derived optical 
depths which can occur as a result of variations in the shapes and sizes of the dust part icle grains (Tab le [TUl) , 
usin g profiles from a var i ety of literature sources, as noted in that table. The reviews bv lHenning et al. I (Epos) 



and iMolster fc Kempeij (|2005l) emphasize that variations in grain size, shape, and metallicity can have an 



i mpact on the shape and peak-wavelength of the silicate profiles. This is well-illustrated in Figures 6 and 7 



of 



Chiar fc TielensI (120061 ) for variations in the dust particles' shapes, including solid versus porous particles, 
and solid glass versus a continuous distribution of ellipsoids. These differences are most pronounced in the 
18 /im region. 

We first explore the impact of the variations in grain shape by fitting three separate amorphous olivine 
profiles to our PKS 1830-211 absorption features: one which is typical of solid spheres, one which is typical 
of a porous continuous distribution of ellipsoids, and one which is intermediate between these two extremes. 
Figure [19] illustrates the visible differences in the fits which can result from these variations. As detailed 
in Table 1111 the best fit results from the profile representing the continuous distribution of ellipsoids with 
a porous structure, although none of the profiles adequately fit our substructure, as discussed elsewhere, 
and all have relatively poor values. Furthermore, as detailed in Table [TTl while these variations impact 
the quality of the fit, they have relatively little effect on the derived optical depth which only varies from 
0.11< no <0.13. 

We next consider the impact of variations in grain size, considering three simi lar particle distribution s 



with particle sizes of 2 /j,m, 1.5 /xm, and those in the Rayleigh limit, all taken from iDorschner et al.l (|l995r ). 
We find that these size variations can have an even more significant impact both visually (see Figure [TO]) 
and in terms of the (Table [TT]), compared with the explored particle shape variations. The best fit over 
the 10 /im range is produced by the 2.0 p.u\ particles, while over the full fitting range the 1.5 /im particles 
provide a better fit. However, despite the range in relative quality of of the fits, all poor in comparison with 
our crystalline olivine fits, the optical depths are consistently in the 0.11-0.13 range. 



C.3. Variations in Laboratory Crystalline Olivine Templates 

Given our results in Sj3]that the PKS 1830-211 z=0.886 absorber silicate profile is best represented by 
crystalline olivine, we explore in this section the effects of varying the chemical composition (i.e. the Mg/Fe 
ratio), particle shape, and material temperature on the derived fits. To implement these tests we examined 
a total of 19 different varieties of crystalline olivines, as detailed in Table [101 The derived fits using these 
profiles are illustrated in Figure 1201 where it is evident that both particle shape and size variations can 
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significantly alter the fit; temperature has a reduced effect over this region of the spectrum, but would play 
a larger role at longer wavelengths. 

We first consider the variations in composition by exploring six distinct olivine compositions, as well as 
two different hortonolite profiles, since there can be variations between different laboratory samples. These 
profiles include both naturally occurring (terrestrially found) compounds, as well as more synthetic blends, 
as described in the references given in Table [TUl As discussed in ij3.3.2[ we see that as the Mg/Fe ratio is 
varied, the leftmost ('^11 fim) and central (^10 /im) substructure features are increasingly well or poorly fit, 
and despite exploring 7 olivine varieties, we do not find any that perfectly fit our data. However, the wide 
range of fits which we do find, lead us to believe that a ratio of Mg/Fe exists that will fully, simultaneously fit 
these two features. T he best fits are obt ained using the two hortonolite profiles, an intermediate Mg/Fe ratio, 
with the profile from ljager et al. ( IQQSl ) providing slightly superior results. The next best fit is provided by 



forsterite, the Fe-free end of the olivine compositional spectrum. As with the amorphous olivine, we find 
similar fits if we fit over solely the 10 /im range, or if we expand fitting to include the full range. The best 
fit produces an optical depth of tio=0.27±0.05, and the full set of profiles ranges from 0.23< tio <0.28, 
well-within the one-sigma limits. Since these templates span the full range of Mg/Fe ratios, we anticipate 
that the true Mg/Fe ratio will likewise produce an optical depth within this range. 

We next explore the effects of variati ons in particl e shap e using four chemically consistent compounds 



(Mgi.gFeo. 18104) close to forsterite, from lFabian et al.l (|2Q0l[ ). These variations, which include considering 



a naturally occurring particle distribution, as well as spherical particles, and two variations of a continuous 
distribution of ellipsoidal (CDE) particles produce even more visually and statistically significant variations 
in the quality of the fit than the chemical variations discussed above. The best fit is produced from the 
powdered olivine, and it is apparent that the CDE distributions are not good fits. However, again despite 
significant variations in fit quality, the derived optical depth measurements are self-consistent within one- 
sigma across both shape variations and fitting regions. 



Lastly, we consider the effects of variations in temperature. As addressed, in detail, by iKoike et al 



(|2006l ). variations in temperature can induce visible shifts in the profile features, particularly at long wave- 
lengths. To investigate whether these can impact our fits, we have considered both fayalite (the pure Fe 
olivine) and forsterite (the pure Mg olivine) over a range of temperatures from 10-300K. The variations for 
forsterite are minimal in terms of optical depth, although Axr,io=0-9. For fayalite, Atio,/u;z=0.02, which is 
consistent within one-sigma, and the io=0-24, which is somewhat smaller. From these tests, we conclude 
that particle shape, size, and composition may significantly impact the quality (x^) of our derived fits, albeit 
not the optical depth, while in this regime temperature variations are largely insignificant. 



C.4. Expanded Subset of Laboratory Templates for Other Minerals 



Lastly, we explore an expanded set of laboratory-derived non-olivine silicate profiles including 3 amor- 
phous silicate minerals, 5 crystalline pyroxenes, 4 phyllosilicates, 3 crystalline silicate blends, and 4 types 
of SiC (Table [TOl with included references for template origins). These substances span a wide range of 
chemical compositions, some of which have been proposed as poten tially significant contributors to Galac- 
tic dust sources, such as the phyllosilicates (jPorschner et al.lll978l ). and others of which are known to be 
prominent in specific environments, such as SiC which is found in meteorites and carbon stars. This large 
variety of compounds are illustrated in Figure [5T] However, we find that none of the materials fit the silicate 
absorption feature as well as virtually all of the crystalline olivine silicates, although as discussed in some 
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of the materials such as sihca, synthetic enstatite, serpentine, and SiC can fit a subset of the structural 
features. In conclusion, while in principle a combination of several of these minerals might be constructed 
to produce a plausible fit to the substructure, and certainly several of them are candidates to explain the 
~9 yum feature in combination with crystalline olivine, they do not provide a compelling alternative, by 
themselves, to crystalline olivine when explaining the origin of the detected absorption feature. 
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Table 1. PKS 1830-211 Spitzer IRS Observations 



Order 


A-Coverage 






to()s("l) 


SL2 


5.13-7.60 


60 


19 


38 


SLl 


7.46-14.29 


60 


20 


40 


LL2 


13.90-21.27 


120 


10 


40 


LLl 


19.91-39.90 


120 


20 


80 



Note. — Summary of Spitzer IRS observations. In 
the first column we list the spectral order, followed by 
the wavelength coverage (in /^m), the exposure time per 
frame (in seconds), the number of exposures, and the 
total on-target exposure (in minutes) for each nod posi- 
tion in the given spectral order. This total time is double 
texpNexp, because of the two nod positions. 
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Table 2. Template Summary: Observational Profiles 



Profile 






Ret 


Description 


GCS3 


9.60 


7.52-34.94 


1 


Galactic center source GCS-3 (Galactic diffuse ISM) 


Trap 


9.53 


8.06-12.83 


2 


Trapezium region of Orion nebula (dense molecular cloud) 


IC5146 


9.78 


5.54-13.97 


3 


Q21-6 G8.5 Ilia Fe-0.5 field star behind IC5146 (quiescent dense cloud complex) 




9.69 


8.00-13.50 


4 


^ Cephei red supergiant (stellar env.) 


AGB 


10.08 


7.04-29.96 


5 


model profile for O-rich AGB star OH/IR 127.84-0.0 (stellar env.) 


UL06301 


10.20 


8.16-33.55 


6 


ULIRG 06301-7934 (extragalactic) 


UL06301t 


10.20 


8.16-33.55 


6 


ULIRG 06301-7934, masking 11.3 PAH complex (extragalactic) 



Note. — Summary of optical depth profile templates extracted from the literature for observed, astrophysical, sources. In 
the first column we list our abbreviated name, followed by the peak wavelength for the 10 /J.m silicate feature, the range of 
wavelengths spanned by the profile, the literature reference, and finally a description of the profile. All wavelengths are specified 
in /j,m. For broad profil es, with significant sub-structure, the w avelength at which the deep est optical depth is achieve d is listed 
fl) 'Spoon et al.l ||2006l) based on data fr omlChiar & Tielen3 (I2006D : f 2) iBowev & Adamsonl | |20Qlh based on 
Chiar e t al. {20n|); (4) iRoche AitkenI jl984|) : (5) IChiar fc Tie lend j2o'oq) based on data 



as the pea k. REFS 
data from [Forrest^et^al 



(1975); (3^ 



from flCemper et al T T2002i) : & f6) ISpoon eTal.l 120061) 



Table 3. Template Summary: Laboratory Profiles 



Profile 



Ref 



Description 



AMORPHOUS 



AmOliv 


9.76 


7.28-33.87 


1 


amorphous olivine 


AmPyr 


9.16 


7.26-33.81 


1 


amorphous pyroxene 


CRYSTALLINE OLIVINE - Mg2a;Fe2_2rESi04 


Fayalite 


11.43 


2.00-199.44 


2 


synthetic fayalite (x— 0.0) 


Forst 


11.21 


1.67-199.44 


2 


synthetic forsterite (x— 1.0) 


Horton 


11.30 


1.67-669.09 


2 


natural hortonolite (x— 0.55) 


Olivine 


11.21 


2.00-199.44 


2 


natural olivine (x— 0.94) 


OTHER SILICATE MATERIALS 


AvgSerp. 


10.31 


7.80-13.28 


3 


avg. serpentine (amalgam antigoritc, chrysotile, serpentine) 


AvgSil. 


9.19 


7.79-13.31 


3 


avg. silica (amalgam quartz, fumed silica, hydrated silica, opal, micronized am. silica) 


Bronzite 


10.49 


8.34-505.90 


2 


natural orthobronzite (pyroxene) 


Enstat-syn 


9.27 


1.67-199.44 


2 


synthetic clinoenstatite (pyroxene) 


GO-green-aSiC 


11.36 


6.67-31.25 


4 


green a-SiC (GO sample) 



Note. — As in Table[2] but listing optical dept h profil es for laboratory sources. For pro files with mul t iple p eaks over the spanned r egion, 
the deepest peak is adopted as Apeafc- R EFS: f 1) ISpoon ct al. (200^ based on data fro mlFabian et al ] | |2001D ; (2) El ger et al ] ||1998D with 
tabulated data provided by C. Jager; (3) iBowey Sz .■^^^amsoni (i2002b|) based on data from [FerraroTI^S^ : & (4) iFriedemanr^^al .1 ||1981|) 
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Table 4. Optical Depth Fits for PKS 1830-211 z=0.886 Absorber 



Profile TlO X%10 Pit-Kangc Tj^n Xlj^ll Fit-Xr-ange 



OBSERVATIONALLY DERIVED TEMPLATES 



GCS3 


0.12±0.03 


15.59 


8.60-12.50 


0.13±0.04 


10.46 


8.00-19.45 


Trap 


0.11±0.02 


13.30 


8.60-12.50 


0.10±0.02 


14.30 


8.06-12.83 


IC5146 


^■-■-^-0.02 


13.09 


8.60-12.50 


0.11±0.03 


11.32 


8.00-13.97 


^Ccp 


0.12±0.03 


14.07 


8.60-12.50 


0.12±0.03 


12.75 


8.00-13.25 


AGB 


12+003 

^■-^^-0.02 


10.62 


8.60-12.50 


0.12±0.03 


7.94 


8.00-19.45 


UL06301 


0.13±0.03 


10.27 


8.60-12.50 


0.13±0.03 


7.57 


8.16-19.45 


UL06301t 


0.13±0.03 


9.42 


8.60-12.50 


0.13±0.03 


7.20 


8.16-19.45 



AMORPHOUS SILICATE TEMPLATES (LABORATORY) 

AmOliv 0.11+l°l 14.79 8.60-12.50 0.11±0.03 9.60 8.00-19.45 

AmPyr 0.10±0.03 21.60 8.60-12.50 O.iotoog 14.19 8.00-19.45 



CRYSTALLINE OLIVINE SILICATE TEMPLATES (LABORATORY) 



Fayalitc 


0.23±0.04 


6.33 


8.60- 


■12.50 


0.20±0.05 


6.74 


8.00- 


19.45 


Forst 


0.27±0.05 


5.09 


8.60- 


■12.50 


0.25±0.07 


6.04 


8.00- 


19.45 


Horlon 


0.27±0.05 


3.73 


8.60- 


■12.50 


0-26_|,,0(i 


4.78 


8.00- 


19.45 


Olixiiii:' 


U.28 ' 


7,:S2 


S,(i()- 


■12. .^.O 


U.2(i 1 U.UT 


ih'>2 


8.0U- 


1!). 15 



OTHER SILICATE TEMPLATES (LABORATORY) 



AvgSerp. 


0.15±0.03 


14.82 


8.60-12.50 


0.14±0.04 


13.23 


8.00-13.28 


AvgSil. 


O.lOiO.OS 


27.66 


8.60-12.50 


0.09±0.04 


23.97 


8.00-13.31 


Bronzite 


0.13±0.03 


13.55 


8.60-12.50 


0.13±0.04 


9.38 


8.34-19.45 


Enstat-syn 


0.14±0.03 


13.39 


8.60-12.50 


0.14±0.04 


8.92 


8.00-19.45 


GO-green-aSiC 


0.28±0.06 


14.30 


8.60-12.50 


0.28±0.09 


10.83 


8.00-19.45 



Note. — Best fits to PKS 1830-211 absorption spectrum using our optical depth template 
profiles. In the first column we list the profile name, followed by the peak optical depth 
normalization factor (an, termed r here and throughout the paper), the reduced chi-squared, 
and the wavelength range over which the fit was performed. All wavelengths are in /Ltm. 
Columns 2-4 correspond to fitting performed solely over the (z=0.886 absorber rest-frame) 
10 jum feature (8.6-12.5 /xm). Columns 5-7 correspond to fits extended to cover both the 
10 jum and 18 /Ltm features (spanning 8.0-19.45 ^im), if adequate profile data are available; 
otherwise a more limited fitting range is utilized, based on limits derived from the template 
profile. The best fit is produced by hortonolite (crystalline olivine), both when considering 
the narrower fitting range and when fitting over the full profile. All of the crystalline olivine 
templates fit better than the alternative laboratory and observationally derived profiles. 
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Table 5. Two-Mineral Fitting of the 10 /im Feature 



Profilcl+Profilc2 ti T2 xl Fit-A^ongc T2/T1 



HORTONOLITE + SECOND SILICATE 



HORTON+SILICA 


0.25±0.05 


'^■'^^-0.02 


3.37 


8.60- 


■12.50 


0.09 


HORTON+AVG.SERP. 


0.26±0.05 


01 + 0. 03a 


3.72 


8.60- 


■12.50 


0.02 


HORTON+ENSTATITE 


0.24±0.05 




3.44 


8.60- 


■12.50 


0.10 


HORTON+SiC'' 


0.27±0.05 




3.74 


8.60- 


■12.50 


<0.02 



AMORPHOUS OLIVINE + SECOND SILICATE 



HORTON+AM. OLIVINE^ 0.26±0.05 O.Olto o^"" 3.75 8.60-12.50 <0.02 
AM.OLIVINE+SiC 0.07±0.02 0.19±0.06 6.74 8.60-12.50 2.49 



^1-(T limit formally extends below 0.0; lower limit is considered to be effectively 0. 

value is at the lowest boundary explored in the calculation, T2 — 0.005, and is 
consistent with zero. We, thus, discount this fit. 

Note. — Fits resulting from bi-variate profile combinations, i.e., 

exp[— (a^Ti^TioTm -f bnr2^norm)]- Li thc first columu, we list the two combined 
minerals, followed by the peak optical depth for the first and second minerals, the 
reduced chi-squared of the fit, the wavelength range over which the fit was performed, 
and the ratio of the derived optical depths for mineral 2 and mineral 1; the nomenclature 
is as detailed in Table [4] All wavelengths are in fim. Hortonolite, which produced the 
lowest Xr Table is used for the first profile in every combination but the last. The 
best fit is the hortonolite -|- silica (;^^^3.37) although it is only a marginal improvement 
over the pure hortonolite fit (x^— 3.73). 
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Table 6. Amorphous Olivine Profiles with PKS 1830-211 Masking 



Profile no Xr.lO Fit-Arangc AA T/„ii Xr./ui! Fit-A^ange AX^hift 



LABORATORY - AMORPHOUS - AmOliv 



fiducial 


11 + 003 


14.79 


8.60-12.50 




0.11±0.03 


9.60 


8.00-19.45 




MASKA 


0.16±0.03 


12.79 


8.60-12.50 




0.14±0.04 


8.36 


8.00-19.45 




MASKB 


0.14±0.03 


6.74 


8.60-12.50 




0.13±0.04 


5.97 


8.00-19.45 




SHIFTMASKA 


0.24±0.04 


2.56 


8.60-12.50 


0.55 


0.19±0.05 


5.00 


8.00-19.45 


0.55 


SHIFTMASKB 


0.21±0.04 


1.41 


8.60-12.50 


0.48 


0.16±0.05 


4.10 


8.00-19.45 


0.45 



LABORATORY - AMORPHOUS - AmOlivGPC 



fiducial 


0.13±0.03 


10.46 


8.60-12.50 




0.13±0.03 


7.19 


8.02-19.45 




MASKA 


0.20±0.03 


4.55 


8.60-12.50 




0.16±0.04 


5.53 


8.02-19.45 




MASKB 


0.17±0.03 


2.39 


8.60-12.50 




0.15±0.04 


4.39 


8.02-19.45 




SHIFTMASKA 


0.22±0.04 


2.78 


8.60-12.50 


0.23 


0.18±0.04 


4.64 


8.27-19.45 


0.25 


SHIFTMASKB 


0.19±0.04 


1.83 


8.60-12.50 


0.15 




3.98 


8.27-19.45 


0.25 



LABORATORY - AMORPHOUS - AmOlivGS 



fiducial 


0.13±0.03 


14.29 


8.60-12.50 




0.13±0.04 


8.84 


8.12-19.45 




MASKA 


0.20±0.04 


8.85 


8.60-12.50 




0.17±0.05 


6.46 


8.12-19.45 




MASKB 


0.17±0.03 


4.07 


8.60-12.50 




16+005 


4.47 


8.12-19.45 




SHIFTMASKA 


f, 97+0.05 
■J-^' -0.04 


2.77 


8.60-12.50 


0.38 


0.20±0.05 


4.50 


8.49-19.45 


0.38 


SHIFTMASKB 


0.23±0.04 


1.27 


8.60-12.50 


0.28 


0.18±0.05 


3.45 


8.34-19.45 


0.23 



LABORATORY - AMORPHOUS - AMOLIV1.5 



fiducial 


l-! + 03 


9.85 


8.60-12.50 




0.12±0.03 


6.81 


8.36-19.45 




MASKA 


0.18±0.03 


5.03 


8.60-12.50 




0.15±0.04 


5.38 


8.36-19.45 




MASKB 


0.16±0.03 


3.15 


8.60-12.50 




0.13±0.04 


4.23 


8.36-19.45 




SHIFTMASKA 


0.20±0.03 


4.00 


8.60-12.50 


0.20 


0.16±0.04 


4.89 


8.61-19.45 


0.25 


SHIFTMASKB 


0.17±0.03 


2.94 


8.60-12.50 


0.10 


0.14±0.04 


4.05 


8.50-19.45 


0.15 



Note. — Optical depth profile fits for four amorphous olivine templates (for template details, see Tables [s] 
and IIQI I. listing the fiducial fits (from Tabic [4l . as well as the fits derived when employing masking of presumed 
absorption/emission features in the PKS 1830-211 spectrum, and simultaneous feature masking and wavelength- 
shifting of the template profile, in order to obtain the best fit. Two broad fitting regions are tabulated: the 10 fim 
region (8.6-12.5 ^m) and the full fitting region (8.0-19.45 /J-m), with the fitting region adjusted as appropriate to 
ensure template coverage. The nomenclature is as in Table with the addition of the requisite shift applied in 
the fifth and last columns for the shift -|-mask fits, and all wavelengths in /j.m. The masks are applied to cover 
the following spectral regions: A (9.3-10.05 ^m; 10.35-10.85 fim) and B (9.3-10.05 /xm; 10.35-11.3 fim). The 
formally best fits result when using mask B and wavelength-shifting of the templates, but these fits are also weakly 
constrained in the 10 /im region and assume that plausible physical mechanisms are present which can explain the 
requisite emission/absorption features. The lowest overall Xr '^^ obtained with the AmOlivGS template, employing 
a shift of 0.23-0.28 /^m, and maximal masking. 



Table 7. PKS 1830-211 Spectral Normalization/Processing Variations 



Name Type Order Norm. H no H xl.io B(10 /xm) B no B xl,io B(full) B T/„ii B xl, 



Adopted Normalization 

Fiducial Chebyshev poly. 3 A 0.27±0.05 3.73 Horton 0.27±0.05 3.73 Horton 0.26t° °g 4.78 



Clipping Variations 

Clip v2 Chebyshev poly. 3 A 0.33±0.05 4.58 Horton 0.33±0.05 4.58 Horton 0.32±0.07 4.45 

Clip v3 Chebyshev poly. 3 A 0.30+° °^ 4.71 Horton 0.30j:° ™ 4.71 Horton 0.29±0.07 4.90 

Clip v4 Chebyshev poly. 3 A 0.32±0.06 5.25 Horton 0.32±d.06 5.25 Horton 0.30±0.07 4.94 



Binning Variations 

Binned Chebyshev poly. 3 A 0.29±0.05 2.48 Horton 0.29±0.05 2.48 Forst 0.27±0.07 5.38 



Normalization Polynomial Shape Variations 



Chebyshev4 


Chebyshev poly. 


4 


A 


0.32±0.05 


3.70 


Horton 


0.32±0.05 


3.70 


Horton 


0.31±0.07 


4.66 


Chebyshev5 


Chebyshev poly. 


5 


A 




8.06 


Fayalite 


0.40±0.04 


6.31 


Horton 


0.41±0.07 


6.98 


SplineS 


cubic-spline 




A 


0.31±0.05 


3.64 


Horton 


0.31±0.05 


3.64 


Horton 


0.31±0.07 


4.73 


Legendre3 


Legendre poly. 


3 


A 


0.31±0.05 


4.01 


Horton 


0.31±0.05 


4.01 


Horton 


0.30±0.07 


4.44 


Linear 


linear 


2 




0.26±0.05 


4.37 


HortonB 


0.25±0.05 


3.86 


FayalT300 


0.51_o,o6 


20.97 


Power Law 


power-law 






0.29±0.05 


5.36 


HortonB 


0.29±0.05 


4.00 


AmOliv2.0 


0.28±0.03 


42.20 



Third-Order Chebyshev NormaUzation Shape Variations 



Chebyshev3 (11) 


Chebyshev poly. 


3 




0.26±0.a5 


3.65 


Horton 


0.26±0.05 


3.65 


FayalT300 


0.34±0.08 


5.65 


Chebyshev3 (12) 


Chebyshev poly. 


3 


A,B,C,D 


0.42±0.05 


7.76 


Fayalite 


0.37±0.04 


6.22 


Horton 


0.40±0.07 


5.70 


ChebyshevS (a) 


Chebyshev poly. 


3 


C,D 


0.37±0.05 


5.22 


Horton 


0.37±0.05 


5.22 


Horton 


0.37±0.07 


4.59 


Chebyshev3 (14) 


Chebyshev poly. 


3 


B,C,D 


0.42±0.05 


7.19 


Fayalite 


0.37±0.04 


5.91 


Fayalite 


n ,,; + 0.06 
U.OO_Q Q5 


5.13 


Chebyshev3 {(5)^ 


Chebyshev poly. 


3 


A,E 


0.20±0.05 


4.38 


Horton 


0.20±0.05 


4.38 


FayalTlO 


0.24±0.08 


5.64 



Note. — Comparison of the explored binning, clipping, and normalization variations applied to the PKS 1830-211 spectrum. Wc list in col. 1 the fitting variation 
name, in col. 2 the shape of the fit applied, in col. 3 the order of the fitted polynomial, in col. 4 the code for the regions through which the normalization was required 
to pass (sec below), in cols. 5-6 wc list the peak optical depth and reduced chi-squared using the given normalization and the hortonolitc template, in cols. 7-9 wc list 
the best fitting optical depth template over the 10 fim region and the associated peak optical depth and reduced chi-squared, and in cols. 10-12 we list the same for the 
expanded fitting range. When fitting the quasar continuum profile, the standard applied masking excluded the 10 fira silicate absorption feature, the 18 /j.m silicate 
absorption feature, the '^6.0 /im H2O feature, and the ^6.8 fim feature, as they appear in the z— 0.886 absorber rest-frame. For the primary analysis in this paper the 
adopted combination consists of un-binned, maximally-clipped data with a third-order Chebyshev normalization with fitting forced to pass (sec note-A below) through 
data points near 19 /^m, and standard masking. The quasar continuum fit was required to pass through specific regions in the z— 0.886 absorber rest-frame, denoted 
in col. 4 as follows: (A) points at the longest wavelengths 18.5 /^m; (B) points between the 10 and 18 /xm absorption features; (C) points just shortward of the 
10 fj.m absorption feature; (D) points just shortward of the 6 fj.m H2O absorption feature; & (E) points in the sub-structure peak near 9.6 fj.m. 



Unconstrained within plausible r values explored. 

The ^^6.0 and 6.8 fj-m absorption features were not explicitly masked in the quasar-continuum fit. 
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Table 8. Expanded Template Summary: Observational Profiles 



Profile 



Rcf 



Description 



EXTRAGALACTIC SOURCES: ULIRGS 



ULoosg?" 


9.80 


7.85-30.19 


1 


ULIRG 00397-1312 


ULoesoi" 


10.20 


8.16-33.55 


1 


ULIRG 06301-7934 


U18443-'' 


9.97 


8.38-34.17 


1 


ULIRG 184434-7433 


GALACTIC DIFFUSE ISM 


GCS3 


9.60 


7.52-34.94 


2 


Galactic center source GCS-3 


WR118 


9.81 


7.10-27.61 


3 


diffuse ISM ilium by WC Wolf-Rayet star WR118 


DENSE MOLECULAR CLOUD/CLOUD COMPLEX 


Trap 


9.53 


8.06-12.83 


4 


Trapezium region of Orion nebula 


Taul6 


9.55 


8.07-12.82 


5 


Taurus molecular cloud ilium, by field star Ellas 16 


IC5146 


9.78 


5.54-13.97 


6 


Q21-6 G8.5 Ilia Fe-0.5 field star behind IC5146 


pOph 


9.81 


7.51-12.29 


7 


p Ophiuci ilium, by SSTc2d.J163346. 2-242753 


T-TAURI STAR EMBEDDED IN MOLECULAR CLOUD 


emb-Tau 


9.37 


8.08-12.83 


8 


T-Tauri star Ellas 7N embedded in Taurus mol. cloud complex 


emb-pOph 


9.47 


8.35-13.29 


9 


T-Tauri stsr Elias 29 embedded in disrupted molecular cloud pOph 


STELLAR MATERIAL 




9.69 


8.00-13.50 


10 


^ Cephei red supergiant 


AGB 


10.08 


7.04-29.96 


11 


model profile for O-rich AGB star OH/IR 127.8-1-0.0 


COMETS 


Bradficld 


10.09 


7.67-13.04 


12 


Bradfield 1987 XXIX comet (r=0.99AU), norm, by 335K blackbody continuum 


Levy 


9.87 


7.68-13.43 


13 


Levy 1990 XX comet, norm, by 266K blackbody continuum 



Note. — As in Table[2] but listing the extended set of optical depth profile templates based on observed, astrophysical sources. 
These two comets were selected from a larger sample of 9 such objects, because they p rov ided the best fits to our s ystem. For ease 
in co mparison, the profiles in Tablc[2]arc repeate d here . REFS: (1) Spoon ct al. ( 200^; (2) ISpoon et a h" (2006} based o n dat a 
fro mlChiar T iclcns (2006); (3) Chiar & Ticlcni 
(5) Bowej^^^Adamso 

(8) Bpwev et al., ( 2003) based on da ta fron i3owpy & Adamson (2001); (9) B owc y ct al . (20 ^0^ ) based o n dat a from Hanne^^^^ 
1995h ; (10) iRochelTAitkenl i 19841): fll) IChiar fcTic lcns (2006) based on data from lKemper et alj i2Q02D; (12)|Hanner et al 



;ni i2006D: ( 4) Bowcy & Adamson ( 200r) based on dat a from Forrest ct aT] i 1975 ): 
from Bowcy ct al. (1998); (6) Chiar ct ar ^ 201l|) : (7) Wa.n Breemen et all i 2011 ): 



I 19941) based on data from lLvnch et"al] jl98g|) : &: fl3) Banner ct al. ( 1991) based on data from lLynch et aU 1 11993) . 

^A nearly identical profile, but with the 11.3 FAB complex masked, has also been utilized in fitting; these additional profiles 
have a "t" appended to the profile name. 
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Table 9. Expanded Observational Optical Depth Fits for PKS 1830-211 Absorber 



Profile 



Fit-A^ 



EXTRAGALACTIC SOURCES: ULIRGS 



UL00397 


0.11±0.02 


12.00 


8.60- 


■12.50 


0.11±0.03 


9.02 


8.00-19.45 


UL00397t 


0.11±0.02 


11.53 


8.60- 


-12.50 


0.11±0.03 


8.81 


8.00-19.45 


UL06301 


0.13±0.03 


10.27 


8.60- 


-12.50 


0.13±0.03 


7.57 


8.16-19.45 


ULoeaoit 


0.13±0.03 


9.42 


8.60- 


-12.50 


0.13±0.03 


7.20 


8.16-19.45 


U18443 


0.13±0.03 


11.73 


8.60- 


-12.50 


0.14±0.04 


8.15 


8.38-19.45 


U18443t 


0.13±0.03 


11.45 


8.60- 


-12.50 


0.14±0.04 


8.01 


8.38-19.45 



DIFFUSE ISM 



GCS3 
WR118 



0.12±0.03 
0.12±0.03 



15.59 
14.92 



8.60-12.50 
8.60-12.50 



0.13±0.04 
0.13±0.04 



10.46 
9.67 



8.00-19.45 
8.00-19.45 



DENSE MOLECULAR CLOUD/CLOUD COMPLEX 



Trap 


0.11±0.02 


13.30 


8.60- 


-12.50 


0.10±0.02 


14.30 


8.06-12.83 


Taul6 


0.11±0.02 


13.62 


8.60- 


-12.50 


0.10±0.02 


14.40 


8.07-12.82 


IC5146 




^■-^^-0.02 


13.09 


8.60- 


12.50 


0.11±0.03 


11.32 


8.00-13.97 


pOph 




^^•-^-^-0.02 


16.28 


8.60- 


-12.29 


0.10±0.03 


16.69 


8.00-12.29 





T-TAURI STAR EMBEDDED IN MOLECULAR 


CLOUD 




cmb-Tau 
cmb-pOph 


0.11±0.02 
0.12±0.03 


13.35 
14.49 


8.60-12.50 
8.60-12.50 


0.12±0.03 


14.30 
13.32 


8.08-12.83 
8.35-13.29 


STELLAR MATERIAL 


/iCep 
AGB 


0.12±0.03 

n -19 + 0.03 
^■-^^-0.02 


14.07 
10.62 


8.60-12.50 
8.60-12.50 


0.12±0.03 
0.12±0.03 


12.75 
7.94 


8.00-13.25 
8.00-19.45 


COMETS 


Bradficld 
Levy 


0.16±0.03 
0.16±0.03 


8.36 
9.94 


8.60-12.50 
8.60-12.50 


0.16±0.04 
0.16±0.04 


8.71 
9.02 


8.00-13.04 
8.00-13.43 



Note. — Identical to Table but for the expanded set of observationally-based optical 
depth profile templates. For ease in comparison, the results for observational profiles in 
Tablc[^arc reproduced here. The best fit is produced by the Bradfield comet in the 10 fim 
fitting range, and by the AGB star and ULIRGs when considering the full fitting range. 
None of these fits arc as good as those found for the crystalline olivine silicates. 



Table 10. Expanded Template Summary: Laboratory Profiles 



Profile ^peak Grange Ref Description 



AMORPHOUS OLIVINE (particle shape var.) 
AmOliv 9.76 7.28-33.87 1 amorphous olivine 

AmOlivGPC 10.06 8.02-29.50 2 amorphous olivine-porous, CDE particle distribution 



AmOliv-GS 


9.85 


8.12-29.74 2 


amorphous olivine-solid, spherical particle distribution 


AMORPHOUS OLIVINE (particle size var.) 


AmOlivl.5 


10.15 


8.36-39.41 3 


amorphous olivinc-1.5 /^m particle size 


AmOliv2.0 


10.51 


8.17-39.46 3 


amorphous olivine-2 pm particle size 


AmOlivRay 


9.77 


8.16-39.73 3 


amorphous olivine-Rayleigh limit particle size 


CRYSTALLINE OLIVINE (Mg2icFe2-2»Si04) - composition variations 


Fayalite 


11.43 


2.00-199.44 4 


synthetic fayalite (x=0.0) 


Forst 


11.21 


1.67-199.44 4 


synthetic forsterite (x=1.0) 


Horton 


11.30 


1.67-669.09 4 


natural hortonolite(x=0.55) 


Olivine 


11.21 


2.00-199.44 4 


natural olivine (x=0.94) 


HortonB 


11.31 


7.80-13.30 5 


hortonolite(x— 0.55) 


Chrysolite 


11.15 


7.84-13.30 5 


oli vine-chrysolite 


SanCarlTlOO 


11.18 


8.00-99.72 6 


San Carlos obs. olivine (T=100K) 


CRYSTALLINE OLIVINE (Mgi.9Feo.iSi04)- particle shape variations 


PowderOliv. 


11.21 


2.00-199.44 7 


ground natural olivine crystals, powder in KBr & PE, sizes <ljLtm 


OlivSphere 


10.67 


2.00-10000.00 7 


spherical particles in Raylcigh limit 


OlivCDEl 


11.19 


2.00-10000.00 7 


CDE(l) particle sizes in Rayleigh limit 


01ivCDE2 


10.92 


2.00-10000.00 7 


CDE(2) particle sizes in Rayleigh limit 


CRYSTALLINE FORSTERITE (Mg2Si04) - temperature variations 


ForstT300 


11.19 


8.00-99.72 6 


,Jona forsloritc T=300K 


ForstT200 


11.19 


8.00-99.72 6 


,Jona forsterite T=200K 


ForstTlOO 


11.18 


8.00-99.72 6 


,lona Ibrstorilc T=100K 




11. IS 


.S. (1(1-!)'). 72 (i 


.Jc'ua [oiXcritc T-KIK 


CRYSTALLINE FAYALITE (Fe2Si04 ) - temperature variations 


FayalT300 


11.40 


8.00-99.72 6 


Jena fayalite T=300K 


FayalT200 


11.38 


8.00-99.72 6 


Jena fayalite T=200K 


FayalTlOO 


11.37 


8.00-99.72 6 


Jena fayalite T=100K 



Table 10— Continued 



Profile 






Ref 


Description 


FayalTlO 


11.37 


8.00-99.72 


6 


Jena fayalite T=10K 


OTHER AMORPHOUS SILICATES 


AmPyr-GPC 


9.50 


7.00-29.52 


2 


amorphous pyroxene-porous, CDE particle distribution 


AmPyr 


9.16 


7.26-33.81 


1 


amorphous pyroxene 


AmSilica 


9.04 


7.17-31.25 


8 


amorphous Si02 


CRYSTALLINE PYROXENE (Mg^Fci-j^SiOs) 


Bronzitc 


10.49 


8.34-505.90 


4 


natural orthobronzite (x— 0.88) 


Enstat-nat 


9.80 


1.67-669.09 


4 


natural orthoenstatite (x— 0.96) 


Enstat-syn 


9.27 


1.67-199.44 


4 


synthetic clinoenstatite (x— 1.0) 


Hypcrst 


10.57 


1.67-199.44 


4 


natural orthohypersthene (x— 0.65) 


AvgPryrox. 


9.33 


7.78-13.32 


5 


avg. pyroxene (amalgam enstatite, bronzitc, hypcrsthene, anthophyllite, diopside, pyroxene-pigconite, wollastonite, augite) 


OTHER SILICATES 


Chamositc 


10.00 


7.81-12.50 


9 


chamosite (phyllosilicate Fe- ■4Al[Si3A10io] [OH]6-nH20) 


Chlorite 


10.10 


8.13-12.50 


9 


chlorite (phyllosilicate (Mg,Fe)6_p(Al,Fe)2pSi4_pOio[OH]8) 


Montmor. 


9.52 


7.69-12.50 


9 


montmorillonite (phyllosiUcate m{Mg3 [Si40io] [OH]2}--p}(Al,Fo- ■ ■)2 [Si40io] [OH]2} ■nH20) 


Serp. 


10.42 


8.13-11.91 


9 


serpentine (phyllosilicate Mge [Si40io] [OH]s) 


AvgSerp. 


10.31 


7.80-13.28 


5 


avg. serpentine (amalgam antigorite, chrysotile, serpentine) 


21Fer 


9.15 


7.77-13.28 


5 


21Fcrraro mixture (amalgam 21 minerals incl. olivines, pyroxenes, serpentines, clays, talcs, silicas) 


AvgSil. 


9.19 


7.79-13.31 


5 


avg. silica (amalgam quartz, fumed silica, hydrated silica, opal, micronized am. silica) 


SiC 


BO-black-aSiC 


11.36 


6.67-31.25 


10 


black a-SiC (BO sample) 


GO-green-aSiC 


11.36 


6.67-31.25 


10 


green a-SiC (GO sample) 


Lonza-aSiC 


11.86 


8.01-16.00 


11 


high temp a-SiC product: Lonza UF-15 sample 


AJ-aSiC 


11.95 


8.01-15.99 


11 


high temp a-SiC product: AJ-a sample 



Note. — As i n Table [3j but listing an extended set of optical depth profile tem plates based on laborat ory sources. For ease in comparison, the pr ofile s in Table [3] are r e peate d 
here. REFS: (l) JSpoon et al.l ^WO^ based on data from lFabian et al.l ll200lh: (2HC hiar fe Tielenj ll2006ll based on dat a fromlHenning et a .1 lll999|l: (31 lporschne~et al.l lll995l);_ 



(4) | jager et al.l 1 1998^ w i thta bulat ed data provided by C . Jage r: (5)lBovyev fe A damson I 2002b) bas ed on data from | Ferrarol Jl98j) : (6) lKoike et aU (|2006|) : (7) lFabianet al.l 
(I2OO1I) : fSl lFabian et aD 1I2OOOI I: (91 iDorschner et al.l lll978t) : f lO") iFriedemann et al.i (.1981 '!: & f 111 iMutschke et al.l lll99gl l. 
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Table 11. Expanded Laboratory Optical Depth Fits for PKS 1830-211 Absorber 



Profile Tio xl,lo Fit-Arange TfuU xljuU FH-Aronge 



AMORPHOUS OLIVINE (particle shape var.) 

AmOliv 0.11+°°| 14.79 8.60-12.50 0.11±0.03 9.60 8.00-19.45 

AmOlivGPC 0.13±0.03 10.46 8.60-12.50 O.lSibO.OS 7.19 8.02-19.45 

AmOliv-GS O.lSibO.OS 14.29 8.60-12.50 0.13±0.04 8.84 8.12-19.45 



AMORPHOUS OLIVINE (particle size var.) 

AmOlivl.5 O.ialool 9.85 8.60-12.50 0.12±0.03 6.81 8.36-19.45 

AmOliv2.0 0.12±0.02 9.35 8.60-12.50 0.11±0.03 7.02 8.17-19.45 

AmOlivRay 0.13±0.03 14.32 8.60-12.50 0.13±0.04 8.81 8.16-19.45 



CRYSTALLINE OLIVINE (Mg2xFe2_2xSi04) - composition variations 



Fayalite 


0.23±0.04 


6.33 


8.60-12.50 


0.20±0.05 


6.74 


8.00-19.45 


Forst 


0.27±0.05 


5.09 


8.60-12.50 


0.25±0.07 


6.04 


8.00-19.45 


Horton 


0.27±0.05 


3.73 


8.60-12.50 


26+" '"' 


4.78 


8.00-19.45 


Olivine 


0.281S:™ 


7.32 


8.60-12.50 


0.26±0.07 


6.52 


8.00-19.45 


HortonB 


0.26±0.05 


4.61 


8.60-12.50 


0.25t°:S« 


5.33 


8.00-13.31 


Chrysolite 


0.28±0.05 


6.04 


8.60-12.50 


0.27±0.06 


6.74 


8.00-13.30 


SanCarlTlOO 


"■^'-0.05 


8.50 


8.60-12.50 




7.56 


8.00-19.45 



CRYSTALLINE OLIVINE (Mgi.9Feo.iSi04)- particle shape variations 



PowderOliv. 


r, r,o-|-0.06 


7.32 


8.60-12.50 


0.27t°:S? 


6.85 


8.00-19.45 


OlivSphere 


19+0°^ 


30.23 


8.60-12.50 


0.221°;- 


18.31 


8.00-19.45 


OlivCDEl 




30.23 


8.60-12.50 




18.31 


8.00-19.45 


01ivCDE2 


0.23±0.05 


14.33 


8.60-12.50 


0.23±0.07 


10.31 


8.00-19.45 



CRYSTALLINE FORSTERITE (Mg2Si04) - temperature variations 



ForstTSOO 
ForstT200 
ForstTlOO 
ForstTIO 


0.27±0.05 
0.27±0.05 
0.27±0.05 
0.27±0.05 


6.52 
7.21 
7.40 
7.62 


8.60-12.50 
8.60-12.50 
8.60-12.50 
8.60-12.50 


0.25±0.07 
0.26±0.07 
0.26±0.07 
0.25±0.07 


7.13 
7.44 
7.61 
7.88 


8.00-19.45 
8.00-19.45 
8.00-19.45 
8.00-19.45 




CRYSTALLINE FAYALITE (Fc2Si04 


) - temperature variations 




FayalTSOO 
FayalT200 
FayalTlOO 
FayalTlO 


0.31±0.06 
0.32±0.06 

(, OO-I-0.07 

g ,o-F0.07 
'^■''''-0.06 


6.17 
6.14 
6.25 
6.38 


8.60-12.50 
8.60-12.50 
8.60-12.50 
8.60-12.50 


29+° °'' 
0.30±0.08 
0.S1±0.08 
0.31±0.08 


5.24 
5.29 
5.33 
5.36 


8,00-19.45 
8,00-19.45 
8.00-19.45 
8.00-19.45 




OTlllili AMOliPliOUS 


SILICATES 






AmPyr-GPC 0.11±0.02 
AmPyr 0.10±0.03 
AinSilica 0,0!) ■ 0,0:! 


14.42 
21.60 

2f<,02 


8.60-12.50 
8.60-12.50 

,S,(j0-12.,^i0 


0.11±0.03 

10+""* 

0,08 ' 


10.40 
14,19 
10,19 


8.00-19.45 
8,00-19,45 

S, 00-19, 1,-, 



CRYSTALLINE PYROXENE (Mg^Fei_»,Si03) 



Bronzite O.lSibO.OS 13.55 8.60-12.50 0.13±0.04 9.38 8.34-19.45 

Enstat-nat O.lSibO.OS 17.86 8.60-12.50 O.IA+IH 11.37 8.00-19.45 

Enstat-syn 0.14ibO.OS 13.39 8.60-12.50 0.14ib0.04 8.92 8.00-19.45 
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Table 11 — Continued 



Profile 


TlO 


2 

Xr.lO 


Fit-Ar„„ge 


Tfull 


2 


Fit-Arange 


Hyperst 


0.12±0.03 


12.79 


8.60-12.50 


0.12±0.04 


9.65 


8.00-19.45 


AvgPryrox. 


0.11±0.02 


13.24 


8.60-12.50 


0.10±0.03 


13.00 


8.00-13.32 


OTHER SILICATES 


Chamositc 


0.14±0.03 


13.71 


8.60-12.50 


0.14±0.03 


12.95 


8.00-12.50 


Chlorite 


0.20±0.05 


15.83 


8.60-12.50 


0.20±0.05 


14.59 


8.13-12.50 


Montmor. 


0.11±0.03 


24.42 


8.60-12.50 


0.10±0.03 


22.44 


8.00-12.50 


Serp. 




17.19 


8.60-11.91 


0.21±0.05 


15.59 


8.13-11.91 


AvgSerp. 


0.15±0.03 


14.82 


8.60-12.50 


0.14±0.04 


13.23 


8.00-13.28 


21For 




"^■^^-0.02 


15.80 


8.60-12.50 


0.10±0.03 


15.02 


8.00-13.28 


AvgSil. 


0.10±0.03 


27.66 


8.60-12.50 


0.09±0.04 


23.97 


8.00-13.31 


SiC 


BO-blaek-aSiC 


0.27±0.06 


13.98 


8.60-12.50 


0.27±0.09 


11.57 


8.00-19.45 


GO-green-aSiC 


0.28±0.06 


14.30 


8.60-12.50 


0.28±0.09 


10.83 


8.00-19.45 


Lonza-aSiC 


0.18±0.06 


27.01 


8.60-12.50 


0.18±0.08 


17.13 


8.01-16.00 


AJ-aSiC 


0.19±0.06 


26.56 


8.60-12.50 


0.18±0.08 


16.95 


8.01-15.99 



Note. — Identical to Table [4l but for the expanded set of laboratory-based optieal depth 
profile templates. For ease in comparison, the results for the laboratory profiles in Table 
are repeated here. The best fit is produced by hortonolite (x^ —3.73 over the 10 fim region 
Xr — 4.78 over the expanded fitting region.) In general the crystalline olivines produce 
the best fits, although slight variations in metallicity/chemical composition, temperature, and 
particle shape can have a marked effect on the profiles and derived fits. 



PKS 1830-211 Spectrum PKS 1830-211 Normalized Spectrum 




Wavelength (/im) Wavelength (/.im) 



Fig. 1. — PKS 1830-211 spectrum. Left: Observed flux versus absorber rest-frame (z=0.886) wavelength 
with the colors denoting individual spectral orders as follows: red-SL2; green-SLl; blue-LL2; and magenta- 
LLl. The black line depicts the 3rd order Chebyshev polynomial fit to the quasar continuum used to 
normalize the spectrum. Right: The resultant normalized spectrum, as a function of absorber- rest-frame 
wavelength, with the dashed line denoting unity. 
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Fig. 2. — Optical depth template profiles for observed (astrophysical) sources drawn from the literature, as 
detailed in Table [21 split into two rows for clarity. The left column depicts the profile over the 10 /im silicate 
feature. The right column depicts the profile over the extended fitting region including the 18 fim silicate 
feature, if data is available. 
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Fig. 3. — Similar to Figure [21 but for the laboratory profiles (Table [3]) drawn from the literature. Note the 
complex structure for the crystalline silicate profiles, in contrast with the the relatively featureless amorphous 
profiles. Small variations in chemical composition can produce notable differences in the profile, as illustrated 
by the 4 crystalline olivine templates in the middle panel. 
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Fig. 4. — Continuum-normalized spectrum of PKS 1830-211, in the rest frame of the z=0.886 absorber, 
overlaid with fits using the observed optical depth template profiles. The left panels show fits performed 
solely over the 10 fim feature (8.6-12.5 /im). The right panels show the fits to the expanded range (8.0- 
19.45 fJ-To), when available; otherwise a more limited fit is depicted. Abbreviations for the source names are 
as in Table [5] None of the astrophysically observed optical depth profiles replicate the observed structure in 
the PKS 1830-211 absorber 10 /im region. 
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Fig. 5. — Similar to Figure IH but for the laboratory optical depth templates. Abbreviations for the source 
names are as in Table |3l The rough agreement between the 10 /im and combined 10+18 /im template fits 
suggests that the continuum normalization constraints applied at long wavelengths are reasonable. This is 
addressed further in Appendix [X] The only profiles which reproduce the strong double-peaked profile within 
the 10 /im silicate feature region are the crystalline olivines. 
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Fig. 6. — Six of the most viable profile fits to the 10 fim absorption feature, overlaid on the PKS 1830-211 
normalized flux curve, in the rest frame of the z=0.886 absorber. The top two rows illustrate the crystalline 
olivine profile fits. Note that while the central profile peak is well-matched by hortonolite, the rightmost 
peak is offset to slightly shorter wavelengths relative to hortonolite. The trend is reversed for forsterite. This 
suggests the need for an olivine composition which is intermediate between the two profiles, or which exhibits 
slight variations in the grain morphology. Formally hortonolite provides a slightly better fit, but none of the 
olivine profiles reproduces the leftmost peak. This lower-wavelength peak could be fit by serpentine (which 
also contributes to the slightly under-fit central peak) or by silica. We, thus, conclude that a mixture of 
olivine, and either a second silicate material or ammonia absorption, as discussed in t j3.3.21 is the most likely 
explanation for the composition of the observed silicate dust. 
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Fig. 7. — Bi-variate fits combining the best-fitting hortonolite profile with other templates, overlaid on top 
of the normalized flux curve for PKS 1830-211, in the rest-frame of the z=0.886 absorber. The best fit is 
produced by hortonolite-f silica, although in every case hortonolite is the dominant component. The minimal 
relative contributions from the secondary constituents produces the visual similarity between the 4 different 
fits in the top panel. We also illustrate the combination of hortonolite-t-amorphous olivine and amorphous 
olivine-|-SiC, for comparison; in both of these cases amorphous olivine provides the more minor contribution 
to the fit. 
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Fig. 8. — Bi-variate fits combining the hortonolite and enstatite (left) and silica (right) optical depth profiles, 
overlaid on the PKS 1830-211 normalized fiux curve, in the rest-frame of the z=0.886 absorber. These two 
combinations formally produce the best fits; although the improvement over the pure hortonolite profile is 
slight. 
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Wavelength (/^m) 

Fig. 9. — Illustration of a test to see whether a version of the PKS 1830-211 spectrum in which the 
substructure features are masked (and attributed to alternative absorption/emission sources) could be well- 
fit by amorphous olivine. We depict these masked regions, which were excluded from the fit, with shaded 
grey areas. The top panel masks two 'emission' features (mask A), while the bottom additionally masks 
an 'absorption' feature (mask B). None of the four amorphous olivine templates adequately matches the 
absorption feature, even with masking, because the amorphous olivine profiles are both too shallow and 
peak at lower wavelengths than our observed broad feature. As in previous figures, the left panel depicts the 
10 /im fit, and the right the full fit; wavelengths are in the rest frame of the z=0.886 absorber. 
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Fig. 10. — Similar to Figure[9l but illustrating the fits when both masking of PKS 1830-211 and wavelength- 
shifting of the template profile are implemented. These fits, particularly those for mask B, are an improve- 
ment over the pure masking. However, while the fit is able to reproduce the underlying broad shape of 
the absorption feature, it requires significant masking and consequently emission/ absorption produced by 
alternative materials, in addition to a shift of '^0.1-0.5 /im, perhaps attributable to grain variations. The 
combination of these phenomena seems unlikely in comparison with the crystalline silicate explanation for 
the substructure. 
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Fig. 11. — Examination of whether the substructure features in the PKS 1830-211 spectrum (depicted in 
both the z=0.886 absorber rest frame and the quasa r rest frame) could originate from atomic /molec ular 
transitions which are found in star-fo rming galaxies ([Smith et al.l 120071 ) , ULIRGs ( Spoon et al. 20061) , or 
active galaxies/ AGN (jHao et al.l 120071 ). We show in grey these lines, taken from ISmith et al.l (|2007l ). with 
the breadth of the feature determined from the unresolved FWHM in the four IRS spectral regions, following 
that paper. While some of the smaller/ weaker features in the PKS 1830-211 spectrum could be attributed 
to the depicted transitions, such as the small p eak in the absorber rest frame at the lo cation of the H2 S(3) 
feature, which is prominent in some ULIRGs dSpoon et al.ll200d: Iffigdon et al.ll2006l) . these transitions in 
combination with broad amorphous silicate absorption cannot explain the broader substructure features. 
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Fig. 12. — Illustration of the impact of variations in the adopted clipping and binning on the quasar- 
continuum normalizations, and resultant normalized flux profiles, for PKS 1830-211, in the z=0.886 absorber 
rest-frame. Left: Comparison of the derived profiles when applying additional or reduced clipping of data 
points, with the adopted clipping shown in the top row. Right: Depiction of the impact of binning (bottom 
panels) versus not binning (top panels). Within each panel, the left plot illustrates the PKS 1830-211 
spectrum, with data point colors as described in Figure [1] The black line depicts our fiducial adopted 
clipping, with no binning, combination, while the red line shows the alternative normalization for comparison. 
The right plot shows the spectrum normalized by this alternative normalization. None of these explored 
variations significantly impact the derived fits. These variations are discussed in Appendix [K\ 
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Fig. 13. — Illustration of explored quasar-continuum normalizations, including higher order Chebyshev 
polynomials, a cubic spline, a 3rd order Legendre polynomial, a linear fit, and a power law fit, and the 
resultant normalized flux profiles, for PKS 1830-211, in the z=0.886 absorber rest-frame, with panels as in 
Figure [T^ None of these explored variations significantly impact the derived fits, as discussed in Appendix El 
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Fig. 14. — Similar to Figure [T3l but illustrating variations in the quasar-continuum normalization shape when 
forcing the 3rd order Chebyshev polynomial to pass through certain regions of the spectrum, as depicted 
by the grey shaded regions. Our fiducial fit employes only forcing passage through points at the longest 
wavelengths, as shown by region A in the figures. While, as illustrated on the right, these variations can 
impact the shape of the normalized quasar spectrum and the depth of the absorption feature, they do not 
qualitatively alter our conclusions about crystalline olivine silicates providing the origin for the feature. 
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Fig. 15. — PKS 18 30-211 spectrum overlaid in grey shading with locations of the PAH features taken from 



Smith et alJ (|2007l ). For featur es below the resolu tion of the given IRS spectral order, the unresolved line 
limit FWHM was adopted from lSmith et alJ (j2007r ). It appears that PAH emission could explain some or all 
of the peak shortward of the 10 fim absorption feature. Likewise, it may be contributing to the rise longward 
of the 10 /im feature, although we note that this is the region also contaminated by the scattered light feature 
discussed in Appendix lB.il We note, however, that the usually prominent 11.3 /im PAH complex is absent 
from our spectrum, or alternatively, that the rightmost absorption substructure valley may be deeper with 
some superposed 11.3 /im emission. 
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Fig. 16. — Comparison of the PKS 1830-211 spectrum with the PAH-dominated, rest-frame, IRS spectra 
of sta r-forming galaxies NGC 5866 and NGC 3198, both noted to be sihcate rich, taken from ISmith et all 
(|2007[). These compari son spectra (shown in blue) have been normalized by the dust-|-stellar continuum 



from lSmith et al.l ([2007) and rescaled/shifted to extend from 0.3-0.8. The black vertical lines demarcate the 
region encompassing the 7.7 and 8.3 /um PAH features, while the light grey vertical lines indicate the region 
spanned by the H2S (3) feature, in the star-forming galaxies. In red, we show these prominent features in 
the star-forming galaxy spectrum rescaled and shifted upward to facilitate comparison with similar features 
in PKS 1830-211. (In the quasar rest-frame panels we omit the red rescaled/shifted H2 S(3) feature.) We 
note the alignment of the 7.7 /im complex PAH feature with the 'emission' feature in the PKS 1830-211 in 
the z=0.886 absorber rest frame, as well as the match in alignment as well as in shape/breadth of the H2 
S(3) feature with the small inflection in our spectrum. We do not see evidence of the 11.3 fim feature in 
our spectrum, or of the lower wavelength PAH features near 6 fim. However, as illustrated by a comparison 
between NGC 5866 and NGC 3198, there can be significant variations in the relative line strengths between 
PAH features in different galaxies. 
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Wavelength (jiir\) 

Fig. 17. — Normalized flux profile for PKS 1830-211, in the z=0.886 absorber-rest-frame, overlaid with an 
amorphous olivine fit to the presumed z=0.19 foreground galaxy silicate absorption. While the purported 
10 /im feature could instead be attributed to the 18 fim feature of the z=0.19 foreground galaxy, the 
corresponding 10 /im silicate feature for the foreground galaxy is absent. Instead, there is an emission- 
type feature, relative to the surrounding absorption, present at the rest wavelength of the z=0.19 galaxy 
10 /im silicate feature. Therefore, we do not believe that the structure we see in the PKS 1830-211 spectrum 
originates from a superposition of 10 /im absorption from the z=0.886 quasar absorber and 18 /im absorption 
from the z=0.19 foreground galaxy. 
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Fig. 18. — Similar to Figure lU but showing the fits for an expanded repertoire of observationally-based 
template profiles, as detailed in Table [H The best fits over the 10 fim feature are provided by the comets, 
ULlRGs, and AGB star, all of which tend to be relatively richer in crystalline silicate materials than the 
other sources. 
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Fig. 19. — Similar to Figure [SJ but showing fits for an expanded set of laboratory amorphous olivine 
templates, as detailed in Table 1101 The top panel shows variations with grain shape, while the bottom 
illustrates variations with grain size. Small variations in these properties can have a marked impact on the 
profile shape, but they are not adequate to explain the substructure dominating the spectrum. 
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Fig. 20. — Similar to Figure[5l but showing fits for an expanded set of laboratory crystalline olivine templates, 
as detailed in Table ITOl In the top panel we illustrate variations in composition (i.e. Mg/Fe ratio), in the 
second panel variations in particle shape, and in the bottom two panels variations in temperature for forsterite 
and fayalite. As illustrated, small variations in composition and particle shape can have a substantial impact 
on the profiles and derived fits. The variations with temperature differences are smaller in this spectral region, 
but become more pronounced at longer wavelengths. 
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Fig. 21. — Similar to Figure 13 but showing fits for an expanded set of laboratory templates, as detailed in 
Table ITOl In the top panel we show amorphous silicates, followed by crystalline pyroxenes, other laboratory 
silicates including phyllosilicates, and finally SiC. In comparison with Figure[20l none of these reproduce the 
substructure as well as crystalline olivine, although some compounds can explain the leftmost substructure 
peak. 



